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ABSTRACT 

The oxidation of acetaldehyde to acetic acid using a 
manganese(ous) acetate catalyst was studied in a single-orifice four- 
inch sparged reactor with a view to modelling the performance of a 
two-phase reactor. 

The kinetic data were obtained in a stirred-tank reactor under 
conditions of high turbulence where chemical reaction was controlling. 
It was found that the over-all rate of oxidation could be described 
by a rate expression which is first order in oxygen, first order in 
acetaldehyde and first order in catalyst concentration. 

In order to model gas-liquid reactors such parameters as gas 
hold-up, interfacial area and physical mass transfer coefficients 
must be known or estimated. A review of the literature concerning 
these parameters is presented along with values determined 
experimentally for the oxygen-acetic acid system. Calderbank's equation 
for interfacial area, Hughmark's equation for gas hold-up and 
Calderbank's correlation for mass transfer coefficients for bubble 
swarms (of average bubble diameter < 0.25 cm) and Hughmark's equation 
for mass transfer coefficients for bubble swarms provide values which 
are in reasonable agreement with the experimental values. 

A simple mathematical model was developed which predicts the 
performance of the sparged reactor satisfactorily under simple operating 
conditions, that is, for the case where the liquid may be assumed to 


be perfectly mixed and a pure gas is used. 
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CHAPTER 1 


INTRODUCTION 

Gas-sparged contacting devices or gas bubble columns, in which a gas 
is bubbled through a deep liquid without mechanical agitation, are common- 
ly used in industrial practice as absorbers, strippers and reactors. Some 


of the current applications of gas-sparged contactor technology are listed 


below (1). 


Process Reactants 

Acetaldehyde Oxygen, ethylene 

Acetic Acid Oxygen, acetaldehyde 
Acrylonitrite Acetylene, hydrogen cyanide 
Ammonium Nitrate Ammonia, nitric acid 
Cyclohexane Hydrogen, benzene 

Ethylene dichloride Ethylene, chlorine 

Hydrogen Peroxide . Air, hydroquinone 
Oxo-processing Hydrogen, carbon monoxide, olefine 
Phenol Air, cumene 

Synthetic hydrocarbons Hydrogen, carbon monoxide 


Fermentations, waste oxidations and the general hydrogenation-oxidation- 
chlorination of hydrocarbons are also carried out in sparged reactors. 
Many reaction situations can occur in sparged contactors and no 
simple set of design procedures are available. A few typical situations 
are the following. 
lis Reactants may be in the liquid phase, with the gas serving 


only as an agitating medium. 
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2a A component of the gas may react directly with a component 
of the liquid. 

3. A component of the gas may react with a component of the 
liquid, the reaction taking place on the surface of a 
suspended solid catalyst. 

4, Two (or more) components of the gas may react with a 
component of the liquid or a suspended solid serving as a 
catalyst but not otherwise participating in the reaction. 

In order to be able to design a two-phase sparged reactor, one has 

to deal with gas dispersion, dynamic gas hold-up, mass transfer and 
chemical kinetics. In the case of simultaneous mass transfer and chemical 
reaction, one must determine the relative contribution of chemical 
reaction rate and mass transfer rate on the overall reaction rate. Rate 
of mass transfer is a function of gas hold-up and average bubble diameter. 
Gas hold-up and average bubble diameter are functions of hydrodynamic 

and geometric parameters. 

A considerable volume of literature is now available that is concerned 
with the subject of two-phase flow. The majority of literature covers the 
area of hydrodynamics and heat transfer and a small amount of work is 
reported in the field of mass transfer with and without chemical reactions. 
Design-information on sparged reactors is scarce. 

The main purpose of this work was to study the performance of a two- 
phase sparged reactor in which a desired product was formed by the 
reaction of gas and liquid reactants and to investigate possible design 


procedures. This involved studies to establish the reaction kinetics 
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determining suitable models for sparged reactors, evaluating procedures 
for the various design parameters and measuring conversions in a four 
inch diameter semibatch sparged reactor. Kinetic data were obtained in 
an intensively-stirred autoclave where an attempt was made to eliminate 
mass-transfer effects by creating a great deal of liquid turbulence. 

The reaction chosen for this study was the oxidation of acetaldehyde 
to acetic acid because of the mild experimental conditions required and 
also because this reaction is industrially important. Chemcell Limited 
manufactures acetic acid by the oxidation of acetaldehyde in a sparged 


reactor and has a considerable interest in this reaction. 
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CHAPTER okt 
OXIDATION OF ACETALDEHYDE 
When exposed to oxygen or air, an aldehyde is slowly oxidized to 
the corresponding peracid, which may further react with the aldehyde to 
form the normal acid. The rate of oxidation is decreased by a long list 
of inhibitors and is increased by light and certain catalysts, for 
example, the acetates of Mn, Co and Cu (2). Further, depending on the 
experimental conditions, e.g. temperature, concentration of liquid and 
gas reactants, type of catalysts, purity of reactants, etc., different 
products are formed. For example, for the oxidation of acetaldehyde, 
the reaction can be controlled to give the desired product as peracetic 
acid, acetic anhydride or acetic acid. In the presence of manganese 
(ous) acetate catalyst, oxidation of acetaldehyde gives acetic acid as 
the main product. The over-all reaction is 


Poe 
CH ,CHO + 20, CH COOH. I AE coe pe et Wee oe Bt 


A. Chemistry 6f-Oxidation .of Acetaldehyde and the Reaction Mechanism: 


Initiation and the Induction period:- 

In the uncatalyzed oxidation there is always an induction period of 
unpredictable length (3,4). Use of a catalyst shortens the induction 
period and in its presence there may or may not be a measurable induction 
period depending on the catalyst concentration. For example, Bawn and 
Williamson (4) observed that with Jone molar Mn (OAc), there was an 
induction period, whereas with 4 x we molar there was not. 

When using Mn salts, the end of the induction period occurs when 
there is a change in color of the reaction mixture from pink to dark 


brown (5). At this time there is a rapid uptake of oxygen and a considerable 
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evolution of heat. The color changes noted are associated with a change 
in the valence of the catalyst from Mn** to Mn" *(6). 
The first reaction in the initiation process is that of molecular 
oxygen with AcH (acetaldehyde) to give an acetyl radical:- 
CU AOR aL yar ax G0 2 s0)5: can a mC POR RRS OEE, 
The acetyl radical reacts further with molecular oxygen to give 
the acetylperoxy radical:- 
CH CO © + 0, 7 CH C000" i Pey PAT ie ee Ty ee ree Oe 
Peracetic acid is then formed from this latter radical either by 
reaction with acetaldehyde (reaction 2.4) or acetic acid (if acetic acid 
is used in the reaction mixture) or by a combination of both. 
CH C000" oe ChCHO 72 CH,,COOOH + CH CO’. ee ee es) 


CH,,C000" + CH COOH) > CH, COO0H t CH,CO0'~ Se at emt ol Ps) 


The HOOAc (peracetic acid) so formed then oxidizes the ionic 
catalyst, e.g. Mn**, to its next highest valence state, either by 
Peactione 2.6) or (2.7). 

CHOCOCOH iM y= CHNCOOs rang 1 .Nn Gnu gis a ae ee C26) 


3 3 


CHUCOUOH Gf Mn > CHeCOC = nO 4M.” 


+ 
3 3 7 


G20) 
The resulting oxidized cation reacts with acetaldehyde (AcH) to form the 
acetyl radical (reaction 2.8) which initiates a chain reaction as 
depicted by equations (2.3) and (2.4) 

CH,CHO + Nn*** + CHACO’ + Mn”) +H. we 1 1 1 ws .(2.8) 
It should be noted that reaction (2.8) is much faster than reaction 
(2.2) and is more important than reaction (2.2) after the induction 


period is over. The reduced state of the cation (as formed in reaction 


(2.8), for example) is re-oxidized via reaction (2.6) or (2.7). 


. J P se 
- 7 = 


From the above reactions, it can be seen that the induction period 
is shortened or may be completely eliminated upon the addition of Mn (OAc ) 
to the reaction mixture (5). Conversely, it is lengthened substantially 
by the addition of benzoquinone. This latter substance is well-known 
for its ability to react with (i.e. destroy) free radicals. Many other 
cations have been proposed for the initiation reactions and the following 
order of activity is suggested (3):- 

Mn(O Ac), > CO( OAc), > Cu(OAc),, 
(1.1x107°M) (6.4x107°M) (50 x 107°M) 

The presence of water in the reaction mixture causes an increase in 
the induction period, even when Mn(OAc )., is added (5). This effect may 
be due to the conversion of the Mn salts to Mn(OOH) which has in fact 
been isolated from oxidation mixtures and found to be inactive as a 


Cabalysts(5)). 


The Chain Reaction and Formation of Peracetic Acid:- 

The following mechanism for the chain reaction is widely accepted 
for both the liquid and vapor phase oxidations (3,4). The acetyl radical 
formed in reaction (2.8) reacts with oxygen as described previously 
(reaction 2.3) and the resulting acetylperoxy radical reacts with AcH to 
give HOOAc (reaction 2.4). Reactions (2.3) and (2.4) together constitute 
the "chain reaction"! - in other words one of the reactants in reaction 
(2.3) is produced in reaction (2.4) and one of the reactants in reaction 
(2.4) is produced in reaction (2.3). The overall effect is) the production 
of HOOAc from AcH and 0,, (reaction 2.9). 


SARA tea ee ee ie eee is, I te ei 
CH CHO + 0, =; CH,COO0H ( ) 
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This reaction will be maintained as long as AcH and 0, are provided, 


2 


and as long as no CH CO" or CH C000" are consumed by routes other than 
reaction (2.3) or (2.4). This latter requirement is not always met 
because of the chain termination reactions (discussed later) and in order 
to sustain the chain reaction, reaction (2.8) is called inte play to 
replenish the supply of acetyl radicals. However, reaction (2.8) can 
function only with a supply of AcH and Mn’ ** ions, the latter being 
obtained by reaction (2.6) or (2.7). 


+++ +++ 
Con "Co ) does not affect’ reaction 


It should be noted that Mn 
(2.3) or (2.4) but is partly responsible for the formation (in reaction 
2.8) of one of the reactants in reaction (2.3) and, as will be shown 
later for the decomposition of HOO0Ac. Therefore, the metal ion 


concentration enters into the overall kinetic rate equation for the 


oxidation. 


The Chain Termination Reaction:- 


CO” ahd 


As mentioned earlier any reaction which consumes the CH, 


CH C000" radicals will inhibit the chain reaction (reactions 2.3 and 2.4). 
The effect of this chain inhibition will be much greater for a batch 
reaction than for a continuous one. Some biacetyl may form as a reaction 
product via reaction (2.10). 
2CH,,CO° = CH.,CO — COCH,. Grea, Vn PRP ROe tare one 5 (25 LO) 
Carbon monoxide and products produced from methyl radicals are 
also noted as oxidation by-products. These would result from reaction 
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It is more likely that the decomposition of the acetylperoxy 
radicals is the main chain-breaking reaction. They can, for example, 
react according to reaction (2.12), 

2CH,CO00° + CH,CO°OO"COCH, + 0,» ee eee ee we (2,12) 


to give an unstable product, which breaks down further to stable 


products likeli CH.4)C0.,)0,., eter  FurthePsrsince a0 inhibits the 


2 2 


uncatalyzed as well as the catalyzed oxidation, it is possible that a 
reaction such as- (2.13) is involved; 
CH,C000° pple le CH ,COO0H ald Olas (O08 weer ey Ach Meee ee sh Dee) 


The “OH radical formed does not participate in the chain reaction. 


Conditions For Maximising Peracetic Acid:- 

Yau et al.(7) studied oxidation of acetaldehyde with a view to 
maximising peracetic acid as the final product and maintained the 
following experimental conditions: 

alls Temperature: 5 - 15.C., so that peracetic acid 

formed in reaction (2.9) is stable. At higher 
temperatures, peracetic acid will decompose 
easily to acetic acid (reaction 2.20, given later). 
ne Catalyst type and concentration: CO(OAC), » 
20 - 40 ppm, which is good for the initiation and 
low concentration is used so that there is no 
further catalytic decomposition of peracetic acid 
by reaction (2.20). 
oe AcH Concentration: 5 - 10% by volume supported by 


industrial practice. 
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Formation of Acetic Acid and Acetic Anhydride :- 

In this section the formation of acetic acid and acetic anhydride 
as the main products from the principal intermediate peracetic acid 
will be discussed. 

(1) The Catalytic Decomposition of Peracetic Acid (H0OAc):- 

It is well known that a solution of HOOAc in HOAc is stable 
for long periods of time (4); however, in the presence of over Lesxio. i 
of Mn(OAc), or over 1.3x107°M of COCOAc), (6) a vigorous reaction occurs 
(accompanied by a color change associated with the appearance of a 
trivalent cation). The products of the reaction are HOAc, 0, and small 
amounts of CO,» CO and gaseous hydrocarbons (5). No acetic anhydride 
(Ac,0) was found as a product (5). It is suggested that the 
decomposition of peracetic acid to acetic acid proceeds by the steps 


C250 apice (7. to](2 19), 


Mae CH,COOOH > Nias ing Oete Upset lit cole gl 25) 
CH,COO” + Game CH COO + in eRe wee ee ee Cony) 
HNC OG mr CHICOOD TE Redmctar tn C2 C000" ca s.ce as eed) 
‘ { i { 5 . e « « 6 . é 2 a 6 
CH,CO00" + CH,COOOH +CH,COD” + 0, + CH,COOH (2.16) 
CH,COOOH + Neel es CH,CO00" + Ve ee re, tena any) 
CH, COOOH + CH,CO00” + 0, + CH,COOH + CH.coo". . . . . .(2.18) 
o 3 2 oO fe) 
CRG es CHECC OS enn Cot ee PO en ee, ace 


3 3 
The overall reaction is given by equation (2.20):- 
5CH_COOOH + Mn’’ > 4CH.COOH + CH,CO0° + 20, 
3 3 3 2 
+ 


.OH SeMec” (2.20) 
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(2ioebirecte Reaction-of Acetaldehyde and Peracetic Acid to Give Acetic Acid 
and/or Acetic Anhydride :- 


The decomposition of HOOAc by AcH was originally believed 
to occur by reaction (2.21) and was studied by observing the action of 
AcH on HOOAc dissolved in HOAc (5). 

CH,CHO + CH,COOOH + 2CH,COOH «6 6 eee ee ee ee 6 (2,21) 

HOAc catalyzes the above reaction. Further studies (3,4) indicated 

that the reaction proceeded in two steps:- 

GL CHO.+ (CH COOH: sPenerider. (GM. 2e4.boun os uC 2te28 

Peroxide + 2CH,COOH. . 6 6 6 ee eee ee ee ee ee 4 (2,23) 
Either one or both reactions (2.22 and 2.23) are catalyzed by acid. It 
has been contemplated (2,3,5) that one peracid and one aldehyde molecule 
react to form a complex (reaction 2.22) which decomposes into two acid 
molecules (reaction 2.23). These two reactions were studied at low 
temperatures (5), in the expectation that while the fast reaction (2.22) 
would be substantially slowed, reaction (2.23) would stop altogether. 
This was indeed found to be the case; and at Shen a bimolecular 
velocity constant of 0.031/mole/em*/sec was found, while at =41°C it was 
0.132/mole/cm*/sec (5). By/calculatien, the velocity constant for 
reaction (2.22) at 20°C was Pel cede en een at least fifteen times that 
of the overall reaction. 

Kagan and Lubarsky (5) isolated a peroxide, m.p. -20°c. Bawn and 

co-workers (4) obtained a peroxide (m.p. 22°c) and found it to have the 
empirical formula C,H, 0 On examining its reactions they came to the 


62 *5. 


conclusion that there was a molecule of "AcH of crystallization" 
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involved, so that the formula became C,H,0,,.CH,CHO . The peroxide was 
decomposed to HOAc solution on the addition of CO(OAc),,. 
Kagan and Lubarsky (5) thought the peroxide had the structure shown 


in reaction (2.24) and that it was decomposed to HOAc according to 


reaction (2.25). 


0 40 Oo” OF 
| er Le 
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\ bg H 
0 OH 
Ge | H 2CH.COOH (2.25) 
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They based their proposed structure for the peroxide on the structure 
of analogous compounds formed by the action of HAO, on AcH or C,H, 00H on 
AcH (4). 

Bawn's group believes that the peroxide breaks down to give Ac,0 
first and that a peroxide of the structure suggested would not give Ac,0 


so readily. They, therefore, proposed a reaction, such as (2.26), for 


the formation of the peroxide. 


orek 0 00H 0 
UY os \ // 
CH,Ce t+ , CCH, + CH,CH - 0 - C - CH,. G26) 
NG iv! 
H \ OH+0 
J 
The product of (2.26) would easily give Ac,0 DyESpALttine (orn H,0 
OOH @) 
, (CH 0 (2. 26a) 
CH.CH - 0-C- CH. - 300) ,0 oe H, aia) ee cate oe 


Carpenter (8) recently studied the catalytic oxidation of AcH to Ac,0 
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and found that the reactions (2.26) and (2.26a) are most favorable at 


fe) : 
56°C and in the presence of cobalt and copper acetate catalysts, 


whereas reaction (2.23) was found to be (93) favorable at 40 - 60°C and 


in the presence of 80 -— 250 ppm Mn (OAc), catalyst. Also water formed 


in reaction (2.26a) must be continuously removed, if ole is to be 


obtained as a major product. This could be done by separating the 


oxidation products into gas and liquid and then by dehydrating the 


liquid products by means of some water-carrying agent; many patents have 


been issued on this process (10). Allen (9) studied the catalytic 


oxidationuer ,AcH to: .AcOH and found that the, main. role of the metal—ion 


(besides initiation) is to catalyze the reaction of Ac 00H with AcH so 


effectively that it becomes the main route to AcOH. The catalytic 


decomposition of Ac 00H is relatively too slow to be of any significance. 


Conditions for Maximising Acetic Acid Production:- 


Ag 


Catalyst: manganese (ous) acetate (2,9,11) 
Temperature: 40 - 60°C (1,9) which increases the rate 
of decomposition of peroxide. 

Catalyst Concentration: 80 - 240 ppm (9) which is 
enough for both initiation and the decomposition of 
peroxide. 

No water in the reactants since water can make the 
catalyst inactive. 

Oxygen content in the reacting gases high: 90% oxygen 
and 10% nitrogen (9,12), supported by industrial practice. 
AcH content low: 2 - 10%, supported by industrial 
practice, probably to lessen the chance of by-products 


formation. 
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Conditions for Maximising Acetic Anhydride :- 


The following conditions should be maintained (10) for obtaining 


a high yield (50 to 80%) of acetic anhydride:- 


fe Oxygen content: 7 - 9% by volume as supported by industrial 
practice. 
oe A mixture of Co and Cu-acetates in certain proportions to 


be established by experiments. 

3. The reaction temperature: about 55°C, SO) thaw tne 
decomposition of peroxide will be faster. 

4. The acetaldehyde content im the redction bath: 4 - 63, 
supported by industrial practice, which will probably 
lessen the possibility of any by-products formation. 

Bis Continuous removal of water formed in reaction 26 (a). 

The Formation of By-products:- 

Some by-products, whose concentration is generally small, may 
also be formed depending on the experimental conditions. 

Carbon monoxide may be formed at 100°C by the decomposition of 
methyl radical, Carbon dioxide is formed at about go°C (13). At low 
temperatures, e.g. a0. = uo°C, the probability of their formation is 
emedl. 13). 

Methanol, formaldehyde and formic acid could be formed (5) 


through the following reaction sequence: 


: LE LC EO he a es ee eel 
CH, ah 0, a CH,0,, ( ) 
aaa Me Rare totes aan. Maen ete Be Se eer ees ort. MN ee 

CH 0, > CHO +00) ans ( ) 
2CH,,0° > CH.,0H ae HCHO Guia? sig hie Rat ease ane. ole apo ah We a eer) 


Formic acid is probably formed by direct oxidation of 


formaldehyde. 
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B. Reaction Kinetics for Oxidation of Acetaldehyde to Acetic Acid:- 


Laboratory stirred-tank reactors have often been used in order to 
obtain kinetics of gas-liquid reactions (11, 14, 15, 16). In this an 
attempt is made to eliminate mass-transfer effects by creating a 
great deal of liquid turbulence and interfacial area (11). Van De Vusse (17) 
has shown that at higher transfer rates the overall reaction rate 
approaches the chemical reaction rate. This is in agreement with the 
concept given by Astarita (15) that when two phenomena occur in series, 
the overall driving force is, in practice used entirely by that 
phenomenon which has a much smaller specific rate. Reaction kinetics 
have to be determined under the conditions when chemical reaction is 
controlling. 

In-aldehyde. oxidation, the order of reaction has been found by 
the method of measuring the initial rate for various initial concentrations 
rather than by studying the change of rate during one complete run 
(2, 3, 11). However, Allen (9) and Yau et al (7) have successfully studied 
the reaction by following the concentration vs time curve. 

Experimental techniques used by other workers, to determine the 
kinetics will be discussed now. 

Since the reaction of acetaldehyde is extremely fast and highly 
exothermic (11), the acetaldehyde is generally diluted with acetic acid 


to give measurable rates of oxidation. 
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Purity of Reactants : 

Aldehyde oxidation is very sensitive to traces of metal ions (2) 
and a large number of other substances and hence it is necessary that 
the apparatus for the purification of aldehyde and the vessel in 


which the reaction is carried out should be scrupulously clean (2). 


The oxygen should be ozone-free and also dust-free since dust 
could be a catalyst for aldehyde oxidation (2). 
Acetaldehyde: 

The acetaldenyde must be purified in the absence of air and 
both the catalysts and inhibitors must be removed from it (2, 4). 
The same source of acetaldehyde, if possible, should be used for all 
the experimental runs so that purity of acetaldehyde is same for 
all the runs (4). 
Acetic Acid: 

Analytic reagent grade giacial acetic acid which should be 
free from peroxides and water content (usually 0.001%) should be 


used. In view of the pronounced effects of water on the oxidation 


fe 


reaction, the same sample or source of acetic acid must be used in al 
experimental runs (4). 
Dissolution Rate: 

Since the reaction takes place between aldehyde and dissolved 
oxygen, it is evident that a quiescent volume of aldehyde may quickly 
exhaust the dissolved oxygen, and the rate of absorption of the 
oxygen from the gas phase becomes controlled by the rate of diffusion 


of the oxygen from the surface into the volume of aldehyde. It is 


quite clear that if the rate of reaction is to be measured, the aldehyde 
must be kept physically saturated with oxygen. The rate of oxidation 
has been shown to be a function of the rate of agitation (2, 16, 18, 19). 
The amount of agitation to maintain saturation depends largely upon the 
volume employed (19), the means of agitation, and the rate of reaction 
as controlled by catalysts and light. 
In the photochemical oxidation of heptaldehyde (2), it was shown 
that a shaking speed of 350 rpm was sufficient to insure saturation of 
2 ml. of aldehyde during photooxidation at a rate of oxygen uptake of 
about 8 ml per min. It is evident that the average surface:volume ratio 
when 2 ml. of aldehyde is shaken is far greater than when 20 ml. is 
stirred. Hence a shaking technique using small volumes of liquid appears 
to insure rapid saturation of aldehyde with oxygen more efficiently and 
conveniently than does rapid stirring of large volumes of aldehyde(2). 
Cooper and Melville (18) conducted experiments to determine the 
effect of varying the stirrer speed on the overall rate of oxygen uptake 
for the oxidation of n-decanal and the results obtained show that the 
oxygen uptake keeps on increasing as the stirring speed is increased 
from 100 rpm to 350 rpm but then above 350 rpm, the uptake is lowered. 
The length of the stirring rod used was 2.5 cm. and its thickness 3 m.m. 
so that the liquid surface was continuously being broken as the rods 
rotated. At the lower speeds, its action amounted to the transmission of 
waves across the aldehyde layer and so the whole sample was being 
subjected to disturbance. At the higher speeds, there was not sufficient 
time during a revolution for the aldehyde swept aside by the stirrer to 


be replaced and so the resultant effect was the formation of a vortex 
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in the film, the surrounding liquid being almost unaffected. 
Unfortunately, the onset of the phenomenon of vortex formation makes 

it difficult to tell whether, at this speed, the observed rates of 
oxidation are really independent of the rate of dissolution of oxygen. 
However, a consideration of the effect of sample size on the rate of 
reaction showed them that this required condition was indeed satisfied. 
It was found that the specific rate of oxygen uptake was independent 
of the quantity of material, other factors, e.g. temperature, pressure 


and stirring speed being constant. 


Selection of Catalyst:- 


According to the patent literature, manganese salts are said to 
be very effective catalysts for the production of acetic acid, whereas 
cobalt salts, a mixture of cobalt and copper salts and vanadium salts 
are stated to oxidise the aldehyde to acetic anhydride in high yield 
(11). The salts of iron, chromium and nickel are less effective 


catalysts. Cobaltous acetate gives high yield of peracetic acid (7). 


Order in Aldehyde :- 


Bowen and Tietz (20) reported that the photooxidation of 
acetaldehyde is first order in aldehyde. However, the rates are 
reported ambigously and it is impossible to tell whether they are 
reported per unit volume of solution or per unit volume of aldehyde. 
The oxidation of benzaldehyde in benzene solution (21) with benzoyl 
peroxide as initiator was found to be first order in aldehyde. The 
cobalt ion-catalyzed oxidation of solutions of acetaldehyde in acetic 


acid (4) was also found to be first order in acetaldehyde. A later 
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paper by Ingles and Melville (22) states that the oxidation of 
benzaldehyde in decane is first order in benzaldehyde, but the authors 
give no data in support of this statement. Venugopalet al (11) 

also found the manganese acetate catalyzed oxidation of acetaldehyde 
to be first order in acetaldehyde. In some exceptional cases (3), the 


order of acetaldehyde was found to be 3/2. 


Order in Oxygen :- 


Almquist and Branch (23) found that the thermal oxidation of 
benzaldehyde is first order in oxygen. The photo-oxidation of 
heptaldehyde (2) was found to be a zero-order reaction between 250 and 
400 m.m and a first order reaction above 400 m.m. At 16°C the thermal 
oxidation of benzaldehyde was zero order in oxygen over a large 
pressure range (16). As the temperature was raised, the reaction 
approached a first-order relation. Venugopal et al (11) found the 


oxidation of acetaldehyde to be a first order reaction in oxygen. 


Order in Catalyst:- 


Bawn and Williamson (4) found the oxidation of acetaldehyde to 
be first order in catalyst concentration, using cobaltous acetate. 
Allen (9) studied the oxidation of acetaldehyde to acetic acid in the 
presence of manganese, cobalt and copper acetates and found the 
reaction to be first order in catalyst. Venugopal et al (11) found 
the oxidation of acetaldehyde to be independent of catalyst 
concentration in the range of 0.2 to 1.6% by weight and they used 


Mn (OAc), as the catalyst. 
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Ethanol changed the order of oxidation of acetaldehyde from 0.5 
to 1 in light intensity (2), anthracene to 0.65, diphenylamine and 
hydroquinone to 0.9. It has been suggested (2) that these differences 
are due to two different mechanisms involved. Adducts changing 
the initiator order to 1 are chain-ending inhibitors; those in whose 
presence the order remains 0.5 are chain-transfer retarders. 
Intermediate orders are due to both mechanisms operating simultan- 
eous ly. 
Rate Equation:- 

As mentioned earlier, the rates of photooxidation of acetaldehyde 
reported by Bowen and Tietz (20) are ambiguous and it is impossible 
to tell whether they are reported per unit volume of solution or that 
of aldehyde. Bawn and Williamson (4) studied the kinetics of oxidation 
of acetaldehyde in presence of cobaltous acetate and hence peracetic 
acid was one of their final products and the rate constant found by 
them cannot be used for oxidation of acetaldehyde to acetic acid in 
the presence of manganese(ous) acetate. Venugopal et al. (11) recently 
studied the kinetics of oxidation of acetaldehyde to acetic acid in 
the presence of manganese(ous) acetate and found the rate expression 


as follows:- 
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However, they found that the reaction rate was independent of 
catalyst concentration when the catalyst concentration was varied from 
O.2 te" loo by weight Pwhich"ic contrary to the results obtained by 
other workers (2, 3, 4, 5). Also they used the volume % of 0, in 
the gas-phase instead of dissolved 0,,~-concentration to determine the 
rate equation and the pressure was not mentioned so that the solubility 
oh 0,, in the liquid mixture cannot be estimated. Further, the values 
of the reaction rate constant given in their data cannot be calculated 
by the rate equation given by them. 


Reaction rate data is scarce and the reliability for design 


purposes of that which is available is questionable. 
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CHAPTER Lil 


SIMULTANEOUS ABSORPTION AND CHEMICAL REACTION: MATHEMATICAL DESCRIPTION 
The differential equation representing the phenomenon of simultaneous 

diffusion and chemical reaction in a liquid phase may be written for any 

reacting gas component C, in the form 

D 3-¢ age 
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in the case of the penetration theory, under the simplifying assumption 


cain 


that the gas-liquid interface is plane,where C = instantaneous concentration 


of reacting gas component in the liquid and r= reaction rate. 

A general solution of the above equation involves a number of 
analytical problems whenever r has other than the simplest relationship 
(i.e. first order) to the concentration of the reactants. However, it is 
possible to find asymptotic solutions of the above equation for very 
simple expressions of r, using the concept of reaction regimes as given 
by Astasita (15). Two time parameters i.e. "Diffusion Time” and "Reaction 
Time" have been defined (15) to assist in recognizing the reaction 


regimes in which a system is operating. 


Diffusion Time:- 


Astarita (15) has defined "diffusion time", th through the equation 
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where , Ko = Physical mass tr. co-efft. 
Ds = Diffusivity of the reacting component into the liquid 
ae diffusion time. 


The physical meaning of the "Diffusion time" is the average life of the 
surface elements, that is to say the interval of time among successive 
mixing processes which make the concentration within the liquid elements 
uniform. The range of values for ty for any practical absorber is 

OFO05. th < 0 Oaeised. Rae. 08 A5ced Fhat. .-. «tla » 0863) 
and for a reasonable value of D, limits of Kk can be obtained from 
equations (3.2) and (3.3) as: 
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Available experimental evidence confirms this range of values (15). 


Reaction Time:- 
The reaction time, tp 


by the reaction to proceed to an appreciable extent and is given (15) by 
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where, C' = Equilibrium concentration of the gas-reactant in the liquid. 


For the case of a first order reaction, 


Depending on the relative values of ty and tp» reaction regimes are 


defined (15) as follows: 
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,» represents the amount of time which is required 
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A. Slow Reaction Regime:- 


The value of the bulk concentration of the liquid-phase 
reactant (b.) is much larger than the value of the bulk-liquid concentration 
of the absorbing gas Cin many practical cases. Except for very fast 
reactions, which will be considered later, it may be assumed that the 
depletion of the liquid-phase reactant near the gas-liquid interface is 
negligible as compared to its bulk concentration. Thus b may be assumed 
equal to bo throughout the liquid-phase, and therefore the reaction rate 
is only a function of C. Let it now be assumed that 

Caen cc! 


th. << th = ee ee ee 0 08.6) 


: . r(C. =O) 

where the reaction time has been calculated at the interface concentration 
level, i.e. the average life of surface elements is much less than the 
time required by the reaction to take place appreciably even if the 
concentration were equal to C. threvghoue the surrace element or the 
liquid. In other words, the reaction may be assumed not to take place at 
all during the diffusion time. 

Thus the equation (3.1) reduces to 
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Astarita (15) concludes that, whatever the kinetics of chemical reaction, 
if the reaction is slow enough for the condition (3.6) to be fulfilled, the 


chemical absorption coefficient K, is equal to the physical absorption 


L 


coefficient ae This does not mean that the reaction does not take place 


at all in the absorber; it does take place in the liquid phase, but it is 


slow enough to be negligible during the short life of the surface elements 
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offiagquid. 


The analysis of such processes is entirely based on the determination 


of the bulk-liquid concentration of dissolved gas. 


Chale 


Let N = average rate of chemical absorption (ML ~T. 


) 


then the absorption rate is given by N = K comics 's Gs and can be 


L 
calculated provided the value of On is known. At the gas-liquid interface, 


the concentration of the absorbing component bees is assumed equal to its 


solubility under the assigned interface partial pressure. 


(a) Diffusional Regime:- 
Reaction rate .r=r (co eC Mb £8) ees? Fon. ease 1) Sov OCS e8) 
If the following condition is fulfilled: 

DRC BAC ee SK EEC! oO). oon cy a, ty eb Gene) 
then the reaction rate per unit interfacial area is much larger than the 
absorption rate. It may be noted that the fulfilment of condition (3.9) 
does not conflict with the fulfilment of condition (3.6) and the 
contemponagy fulfilment of both conditions implies that 
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This condition may easily be fulfilled in any practical absorber (15). 

In such a case, the reaction rate is high enough to keep the unreacted 
gas concentration, Coo practically equal to the equilibrium value C' and 
then in this regime, 

1 al} aie: 
— ame = O 


If condition (3.9) is fulfilled, the rate controlling phenomenon is 
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diffusion, so that the overall driving force is completely used up by the 


diffusion process. This condition has been called the "Diffusional 


regime" by Astarita (15). 


The 


absorption rate in the diffusional regime is given by: 


WS i) m Rtas, 
N Ke a.(C. (ome 


re. 21) 


In this regime, though the chemical reaction does not influence the 


7 gs peers ; : 
value of K. , it keeps the bulk-liquid concentration down to its 


L 


equilibrium value C', although absorption takes place at a finite rate. 
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Kinetic Regime:- 


This regime occurs when the chemical reaction rate is low 


enough so that: 


as 
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The 
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overall driving force is entirely used up by the reaction, 
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oO oO 


coal) 


so that, in practice the liquid phase is everywhere saturated with the 


absorbing gas, C. = C 


is* 


The 


' 
O 


absorption rate under these conditions is given by: 


Nea = v(C '-C') 
fe) 


(8.13) 


This equation shows that the total absorption rate in this regime 


e@) 
(ii) 
(A014) 
(iv) 


(v) 


independent of interfacial area 


proportional to liquid hold-up 
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independent of K, 


proportional to r, and 


influenced by the overall driving force on = ¢', 
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It is important to note that only in the kinetic regime is the rate 
proportional to the liquid hold-up. This fact can be used to determine 
when one is operating in the kinetic regime. Operation in this regime is 
a widely accepted procedure for obtaining kinetic data for gas-liquid 


reactions (15) and has also been used in this work for the same purpose. 


(c) Intermediate Regime: - 
This regime occurs when the chemical reaction rate and diffusion 
rate are of the same order. For a first order reaction and under the 


quasi-stationary hypothesis, the absorption rate is given by 
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KS (mia oo as incorrectly taken by Schafteilen and 
Lee 
KR 


Rueee |) (20), kK is still equal to kK on being in the slow reaction regime, 
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irrespective of the subregimes. 


For any other order n, O<n<], the total absorption rate is given by 
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B. Fast Reaction Regime:- 


In this reaction regime, tp<s<ty and it is assumed that the 


reaction is fast enough to keep the bulk-liquid concentration equal to 


its equilibrium value, C ey 
Vv Ca 
Ne ee OD Ie a G rete Aen, oe en en ena 
G 
O 


(for penetration theory), which for a first order reaction reduces to 
N = We ! 1 = 
N vDjK, (C. Ch) by whieh Ky vD,K, 


For any reaction order n, the absorption rate is given by (15) 
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oF Transition From Slow To Fast Reaction:- 
In this case, tp - ty and the mean absorption rate for a 


first order reaction is given by (25) 
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where t* = total time elapsed from the moment the surface element considered 
has been brought to the surface and S = Danckwert's model parameter, rt, 


bor an N-th order reaction , 


SS ees nau 
Ke (Co oe) + AD eee eee? Sok eG: ae et RR ORS 


Where A is an integration constant and can be calculated as a function of 


- as shown by Astarita (15). 


De Instantaneous Reaction:- 
In this case, the absorbing component and the liquid-phase 
Component Cannot co-exist in the same region of liquid. The condition 


tp<<tp is fulfilled and whenever C # 0, b = O and vice-versa. The 


7% Do 
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i] b) 
tp qc, 


where q = stoichiometric coefficient, is fulfilled. The instantaneous 


absorption rate is given by 
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Et Transition From Fast To Instantaneous Reaction:- 
No analytical solutions are available, though numerical methods 
can be applied under very crude simplifying assumptions (15). 
All the absorption rate expressions above were derived from 
Penetration theory. 
The following table (15) summarizes the conditions necessary for any 


chemical absorption regime: 
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CHAPTER IV 


DESIGN OF TWO-PHASE REACTORS 
The two-phase tank-type reactors of interest can be classified in 
three groups (24) as follows: 
(1) Continuous flow tank reactors (CFTR), in which both the liquid and 
gas flow through the reactor, 
(2) Semiflow batch reactors (SFBR), in which only the gas flows through 
the reactor’and the liquid is held in the reactor. 
(3) Batch reactors, in which both phases are charged in the reactor and 
there is no net flow of either phase from or to the reactor. 
The general mass balance equation for the reacting gas 
component Cin the liquid-phase in any of the well-stirred tank-type 
reactors listed above and operating in the 'Slow reaction regime’ can 


be represented by the following equation:- 


€.8 


fe) ere! 
Qcc, = Cc? ap Kq-Pav, (Y, = eo = PoVy, = ar V1: °5): Teen evil, 1) 


: ‘ : ee Popak 
where a = interfacial area, per unit volume of liquid, L/L. 


Of =-inlet feoncentration of Chin the Mquid Moles /voel. 


Aare 3 
Ve = ‘volune Sof Wliquid phase y-L 


L 
Yo = mole-fraction of component’ C in the gas-phase 
Cl) = bulk=Liquid concentration of CC; Moles/vol. 
° 
Ka = overall mass transfer coefficient, mole/area-time-pressure. 


For a system operating in the 'Fast Reaction Regime’, the reaction 


is fast enough to take place appreciably during the life of surface elements 
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and the concentration oe of the reacting gas-component (C) can be 
approximated to zero (15,27) and there is no accumulation of the gas- 
component C in the bulk liquid. Hence the general mass-balance for 
the gas-phase reactant in the liquid phase can be written as:- 


Q(c. -~ 0) + K -aeP.V CY, 0) ' o ets dO). mection oh.. eres Ca?) 


G C 


where N, = Average chemical absorption rate of the gas-component 
C, Moles / volume-time 
If the gas-component C reacts with the liquid-component b such 


that one mole of C reacts with one mole of b then 
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where x = rate of depletion of b by chemical reaction. 
The general mass balance for the liquid-phase reactant can then 


be written 


A d(b v,) 
Q(b. = bd.) - No: Vi, Fu a eaete Corrs oe a ee Cas) 


The equations (4.2) and (4.3) are also valid in the instantaneous 
reaction regime and the mass flux a can be obtained from the last 
chapter. 

The liquid-phase behavior of any of the two-phase reactors can'be ' 
obtained, in the slow reaction regime, from equation (4.1) and in the 
fast and instantaneous reaction regimes from equation (4.2) or (4.3) by 


deleting the appropriate terms. Fora steady state continuous flow tank 


reactor, equation (4.1) is applicable if the time-derivative is set 


Ou 


equal to zero. For a batch or semiflow batch reactor the terms containing 


the imput and output flow rate,Q, are eliminated. 
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Schaftlein and Russel (24) incorrectly presented equation (4.1) as a 
general mass-balance equation irrespective of the reaction regime. They 
did not seem to realize the difference between the "mass flux across the 
interface" and the "absorption rate due to chemical reaction in the bulk 
phase". 

Several different models have been employed to describe the gas- 
phase behavior. 

In the absence of mechanical agitation the bubbles of gas will rise 
in a swarm through the reactor and their spatial variation over the reactor 
length must be considered. The total mass of the component C in the bubble, 
PV LY /RT is dependent upon both z and t. For a single bubble the 


differential change in the liquid (z, »t,) and(z t,) is equal to the mass 


Oz 


of C transferred to the liquid, “Kaa PV - CH/P), and since the mass 


of C in the bubble is a function of both z and t, 


d(PV, Y/R) a dz ee Dae 
dt Meo on eat ot 
since UL = sy, the following equation results: 
PV PV 
CH 3 b CDP Sa ke 0) 
2 a )t = a =—— “= Cee ae eee 
Kean So — “B) - Uae Se’ Yo) “at Rr a ee 
where ao ratio of surface area to volume of a single bubble 
ie = volume of a bubble 
U, = rise velocity of the bubble 


If there is mechanical agitation and deliberate dispersion of the gas, 
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it is possible that the gas-phase may behave as if it were well-mixed; 


an unsteady state mass-balance for this situation yields a second gas-phase 


model: 
d hie 
Caveoi- GAY. = V - =— a 
ie oto Kaa plPYo CH ) art (N,V. V,) RT ee mere er Reset) 
where , ave = total surface area for mass-transfer 
NOY = total gas volume 
Ka = overall gas-phase mass-transfer co-efficient, moles/ area-time- 
pressure. 


These two equations (4.4) and (4.5) represent the limiting configuration 
for the gas-phase and should provide, when combined with the liquid-phase 
equation of interest, two-phase models which can be solved to indicate 
performance limits for a reactor in which the gas-phase doesn't behave 
"ideally". 

Schaftlein and Russel (24) have summerized the two-phase model 
developments for three reactor types: (1) CFTR, (2) SFBR and (3) Batch 
reactors, and have also discussed the various possible cases for each of 
the reactor types. Only a semi-flow batch reactor (SFBR) and its various 
possible cases will be discussed here. 

Semiflow Batch Reactor: 

In this case liquid is charged in the reactor and gas flows 
through it. Because of this, the liquid concentrations are time-dependent 
but spatially invariant for constant gas flow rate. The gas-phase may be 
described by either the plug-flow or well-mixed models. 

Case A: Plug-flow Gas-Well-Mixed 
Batch Liquid: 


The liquid is assumed to be well-mixed and the general mass 
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balance equation for the gas-phase (4.4) can be written as 


PV. 
CH d b 
V et ae a oe, ie 
Koay pte a, + U, = or : Yo) ee er eae ce ee eee nC te) 


Which can be combined with the mass-balance equation for the liquid phase 
to give overall model equations. Special cases for the SFBR are considered 
below, assuming a first-order reaction. 


Case A-(1) Ve, = constant, first order reaction. 


Here it is assumed that the bubble volume remains practically constant. 


This will be true if the gas-phase reactant concentration, Yoo is small or 


if the amount of mass transferred from the bubble is small. In this case 


then a and UL are also constant and equation (4.6) can be integrated to 


give 


ec3 CH | 
rete OG er ae ban oC ie ZU) ae a oY oe seeds Ga) 


The gas-phase concentration is a function of z and therefore it is 


necessary to define an average value of Y, for the driving force of the 


C 


liquid-phase model, i.e. 
=m of" . One) 
tC akca f Y(zZ) G24 800 Fhe tote yielace se. 


The average values ay and He can be obtained in the same manner. Thus 


Yo is obtained (24) from equations (4.7) and (4.8) as: 


= ie. Vile dae ol 
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Diear ean Beate oteety eget Crt) 


where ny = Ka, RIL/U, 


By substituting equation (4.9) for Ye in the liquid-phase model of 
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interest, the solution for the two-phase model can be obtained. 


Case A-(2) V= constant , Yo = 1, pure gas-phase, first order reaction. 


In this case the bubble volume changes significantly because of 


Mass transferred’ from the bubble into the liquid) "Beth a, and ii are 


functions of the distance the bubble has travelled, z. For this situation 


equation (4.6) can be solved for Vi (2), ie a and U, can be expressed as 


PING LONS. OF V5° 


From the definition of 2 and the assumption that all bubbles are 
spherical, a, can be expressed as follows (24): 
1/3 ee anes 
6 T : 


a, = SY, where Ss. = Cis. er) 


If the rise velocity of bubbles through liquid is assumed to follow 


the Davies-Taylor equation discussed in Davidson and Harrison (26), then 


1/6 
= ' 4.11 
ot ( ) 
1/6 
where W = 0.711 aoe (=) Oe ee ee ae ee ee eam re I oe gD) 


The gas-phase model (equation 4.6) can be integrated after the proper 


substitution (from equations 4.10 to 4.12) to yield V (2) (24): 


Z 1/2 _ ay ppl ~ yy} ? : (i 18) 
ve = NES aK RIG (call poz. LP atmos Pat etter ete : 
a and ee can be obtained as in equation (4.8): 
Pp 
= ie 2 
= err ae aa ae. Vi Inte? ° ° ° ° © . ° e e ° ° 6 ° ° ° (OMe 20) 
2 is ( 3 ob ) 
S] 
CH 
= 4 Pra Sa me EC) ° tn) ° ° « ° ° ° ° © . ° ° © ° © ° ° (4.15) 
id KG RTE Ci 5) 
-1/3 
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These previous results can then be substituted in the liquid-phase 
model of interest and solution for C(t) can then be obtained. 


Case A-(3) V,# constant , Yo 7 1, first-order reaction. 


In this case the change in bubble volume cannot be neglected. The 


overall gas-phase balance is then: 


CH d b 
V. = eS ee (Wy oS 
Koa? pte >) ar eae ( RT Oh hee siete een Pek (ete) es aT Heb) 


Equation (4.17) and equation (4.6) then give: 


CH . 
Ka RT(Y, - =>) + , = Oe ee ee ee ae oe ee ee Ce mI LISE)) 


This gas-phase model can then be solved for YQ(Z) and then average 


values of ai: Vi, and Yo 


to solveimor CG) or b(t). 


must be used in the liquid-phase model of interest 


Case B: Well-mixed Gas-Well-Mixed Batch Laud = 


In this configuration both phases are intimately mixed by 
mechanical agitation and the action of the gas on the liquid. The 
concentrations in both phases are assumed to be time-dependent but 
spatially invariant. 


The gas-phase model equation is then 


CH 
=< ~ i eee = O ° ° rc o ° ) e o ° o © uy, 19 
eee GY Kav, PCY, 5) ( ) 


This gas-phase equation is for the steady state but is coupled with 
the time-dependent liquid phase concentration C(t). 
This equation can be solved for Y,, substituted in the liquid-phase 


equation of interest, and the resulting differential equation solved for 
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Flow Patterns in Two-phase Flow and their Import ance :- 


Two phases flowing together in a tube assume various different flow 
patterns or regimes depending on the relative flow rates, pressure drops 
and tube orientations. The actual conditions which cause the existence 
of individual flow patterns are not completely understood (27). Indeed, 
the transition from one flow pattern to another is difficult to identify. 
It is customary however to identify the flow patterns in the manner 


indicated below (27). 


Flow Patterns 


Horizontal and Helical Increasing Vertical 
Coiled Tubes Gas Flow Tubes 

Bubble Bubble 

Plug Slug 

Stratified Semi-annular 

Wavy Annular 

Slug Dispersed 

Dispersed 


The well-known Baker plot first developed for horizontal low pressure 
two-phase flow (28) and subsequently for vertical and helical coil flow 
data (29) is the most widely used method for mapping flow patterns (27). 
This chart can be used in conjunction with the reactor model to decide which 
flow regime should be most profitably used for particular types of chemical 


reaction. In certain flow patterns, for example, in slug or plug-flow, the 
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hold-up at a given point fluctuates with time. In these cases, time- 
averaged values are used. 

Depending on the flow rates and phase properties of liquid and 
gas, the flow-patterns in the reactor vary and so do the various model 
parameters, e.g. interfacial area, mass-transfer co-efficient, hold-up, 
etc., affecting the reactor performance directly. Details of various 
flow patterns and their effects on two-phase reactor performance can be 


found in the papers by Cichy et al (30) and Roseheart et al (27). 
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CHAPTER V. 


PARAMETER EVALUATION 


In addition to kinetic data one must have values for the physical 
mass transfer coefficient and the interfacial area. The interfacial 
area can be estimated from information on gas hold-up and average bub- 
ble size. This section deals with the bubble behavior and mechanics 
of bubble formation as well as the evaluation of mass transfer 
coefficient and hold-up. The discussion is limited to consideration 
of the bubble regime. 

A. Bubble Behavior and Mechanics of Bubble Formation:- 

Although a single orifice would not normally be used to form 
bubbles in any practical piece of equipment, an understanding of the 
process of bubble formation at a single orifice is a preliminary to the 
study of disperison devices containing multiple orifices. 

The simplest mechanism of bubble formation is that in which the 
bubble is formed very slowly at the open end of a tube immersed vert- 
igally in a digquid... In this: case theybubblée willysgrow untalats 
buoyancy exceeds the surface tension forces tending to hold it onto 
the tube, when it will detach. Assuming the bubble to be spherical, a 
force balance considering only buoyancy and interfacial tension yields 


the following expression for initial bubble diameter d 631): 
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where, qb = initial bubble diameter at the orifice 


d = orifice diameter 


= density of liquid and gas respectively 


= Suntace-tensien of liquid 


g = acceleration due te gravity. 

Datta, Napier and Newitt (32) obtained a large amount of 
experimental data for the air-water system. Working with orifice 
diameters ranging from 0.022 - 0.519 cm (made from glass capillary 
tubes) they confirmed that the above equation was approximately 
correct. They also found that the bubble size from an inclined orifice 
was smaller than a vertical one. This could be explained by the 
tendency of the buoyancy forces to drag the bubble across the orifice 
thus allowing the surface forces to oeprate over only a portion of 
it (33). 

As the gas flow is increased, the bubble diameter at first remains 
| reasonably constant while the frequency of bubble formation increases. 
Then at higher gas flows the frequency of formation becomes approximately 
constant while the bubble diameter increases with gas flow rate. At 
still higher flows bubbles of a single size are no longer produced; 
there is a considerable spread of bubble sizes with production of numerous 
small bubbles (33). Most research workers e.g. Coppock and Micklejohn 
(34), Benzing and Myers (35), Leibson et al. (36), Davidson and Amick 
(37) agree with this behavior of the bubbles with respect to gas flow. 
Bubble Formation Regimes :- 

Based on the gas-flow rate and hence the orifice Reynolds number, 


four bubble formation regimes are defined (33):- 
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(a) Constant Volume Region: - This region has been defined to extend 
up to an orifice Reynolds number of about 200 (33). In this region 
both gas rate and liquid viscosity have negligible effect on bubble 
diameter and equation (5.1) would be expected to hold approximately. 
(b) Slowly Increasing Volume Region:- Davidson and Shuler (38) 
carried out experiments in this region and the following equation is 


suggested for the bubble diameter. 


d= Gonst. 4 


ae {Gu} oR BDAY Ges 


esa) 
where, G = volumetric flow of gas 
and ee viscosity of liquid. 


This is a modification of equation (5.1), which takes into account the 


effects of viscosity and gas-rate on bubble diameter. 

(c) Constant Frequency (Laminar) Region:- As the gas-rate increases 

to correspond to orifice Reynolds number in the range of 1000-2000, the 
frequency becomes independent of gas flow rate and approaches a constant 
value depending on orifice diameter. It ranges from 15/sec , for large 
diameter orifices (35) to 45/sec, for small capillaries (39) being about 
20/sec, for 1/15-1/4 inch diameter vertical orifices (40); the bubble 
volume then increases with increasing gas rate. Only one correlation 
has been developed for do (4) which is claimed to hold for gases at 
orifice Reynolds numbers between 1000 and 2000: 
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Although the observed size of bubbles formed at an orifice is 


reproducible, considerable coalescence and bubble break-up begins to 
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eccur in this region so that equation (5.3) may not give a good estimate 
of the bubble diameter in a bubble column. The above behavior is 
Supported by Davidson and Amick (37), Benzing and Myers (35), Calderbank 
(40), Leibson et al. (36), and Valentin (31). The fraction of the small 
bubbles gradually increases with increasing gas flow rate. 


(d) Turbulent Region: - At gas flow rates corresponding to an orifice 


Reynolds number greater than about 2100, gas bubbles break up due to 
turbulence. This transition point cannot be taken to be a sharp one, 
Since a certain amount of bubble break-up may take place even at a 
lower Reynolds number. Also in this region bubble coalescence begins 
to take place closer and closer to the orifice and a variety of sizes 
of bubbles are formed. Leibson et al. (36) have classified another 
region above a Reynolds number greater than 10,000, where turbulence is 
considered to be fully developed, and the gas jet disintegrates into a 
range of fine and coarse bubbles near the orifice. Siems (42) has 
named this zone as 'jet gassing'. However, according to Rennie and 


Evans (43) even at Reynolds number of 40,000 there is still bubble 


formation followed by bubble break-up, the bubbles now being all toroidal. 


Leibson et al. (36) measured the average bubble size for the air- 
water system by a photographic method and found a sharp fall over the 
Reynolds number region 2100-10,000 from the 'laminar' values (given by 


equation 5.1) to an approximately constant value given by 
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Rennie and Evans (43) have approximately confirmed the above equation. 


Although this equation applies strictly to Single orifices operating at 


10,000, there is reasonable agreement in the range 2000-10000 tor the 


frothing systems which have been studied. 


Bubble Swarms:- 


Koide et al. (44) give the following equation and limits for 


cal- 


culating the average bubble diameter of bubble swarms in air-water system, 


generated from perforated plates 
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WZ 
where N = Weber number = a 
We 6 
2 
N = Froude Number = u 
e Sy 


0 = density 
6 = pore diameter 


w= gas velocity through orifice. 


(5.5) 


The average bubble diameter of bubble swarms from porous plates generally 


ranges from 0.05 to 0.5 cm. and may be estimated by the equation given 


below (44). 


N 
lei) (0) 278 je0) \ as 
= 5 c=) nee 
dibs 3 et ore 


We 


At the present time it appears that the initial bubble diameters 


generated in bubble column are not easily estimated unless one has a 
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fairly low gas rate through the orifice. Experimental studies may have 


i 
i e 


sdoitavps svede edt bemitaes ylesemixonaye ¢ 
tN Tie . 


$e gatterego 2esitizo efgniz ov vi. i 2. & 


ett to2 goo, or=090s egasx efit nt iaBaBsaa8 
.beibusa need ¢ ei dotitw < 


~fe4j sot asimil bas noftsupe gniwollot sat ovis (4H) sis 18 ot 


sMetave Stew-tis ‘of snrmewe efddud to setsmedb atdae SRetsve edt 3 


astslq halaueteet: most B 


u 
e\I-( 98), £002 00M giosgS ‘oe 
Z, | 

pe 


a 
(2 2rd ’ « w : e ‘ss ._ -« « «> ‘s. © . #¢ > ~ a * 
: xT org 


iN 
943 _ ; : ‘oa = a etka 
; sedan vedsW = wv ia’ 


*y = ysdmU shuor 


2% 


wa 
v= 


yvtieasb 


it 
2 
1, 


wstsmeib sveq = 3, 
.soltito dyuowds yrisolev asa = 0 | 


yileianes 2stefq <vo1sg moat Eombve olddvdite | resis if feat 36:18 
avis nettanps: efit Yd betsmites sd yem bas m9 2.0 of 20:0 sort 20 


to be conducted to evaluate mean bubble diameters unless the specific 
system of interest has been investigated and reported in the literature. 
A comprehensive listing of literature sources related to bubble motion 
is presented by Baker and Chao (45), Peebles and Barber (46) and 
Galore (47). 
Bubble Size Distribution:- 

Houghton et al. (48) examined a large number of bubbles in both 
fluidized and foam beds and determined the bubble size distribution. 
The bubble sizes followed a typical probability distribution about the 
most probable value, which was found to vary with plate porosity, gas 
velocity and the properties of the liquid phase. The relationship between 
bubble size and the cumulative percentage of bubbles of a smaller size 
than this, when plotted on log-probability paper is linear, which confirms 
the presence of a/probability function. Observed bubbles smaller than 
the most probable size may have been formed at very small pores, or 
may have been broken off from other bubbles; those larger than the most 
probable size could be produced either by large pores or by coalescence. 
Coalescence seems less probable in the main body of the liquid than close 
to the bubbler plate, since most bubbles seem to coalesce by being pulled 
into the slip stream of the bubble ahead. Coalescence of bubbles moving 
in opposite directions appears unlikely, since the moving films of liquid 
around the front and sides would appear to repel them. Gal-Or (49) and 
Ho and Prince (50) have analyzed the effect of taking bubble size 
distribution into account when calculating the overall physical mass 


transfer coefficient from empirical correlations for single bubble and 
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have concluded that the effect is negligible if the volume surface mean 


diameter is used for calculating the mass transfer coefficient. 


Bis Mass Transfer Coefficient:- 
The usual combined resistance approach is employed to define 
the overall gas-phase mass transfer coefficient, Kea in terms of gas- 


phase and liquid-phase resistances, i.e. 


er on, Tine ee ana Ray ee ae ea 
Ko i Ne 
where, K = individual liquid-phase chemical absorption coefficient, 
length/time 
K.. = individual gas-phase mass transfer coefficient, mole/ 


area-time-pressure. 
In any 'bubble-type' gas-liquid reaction system in which a pure 
gas is employed, the gas-phase resistance, 1/Kg can be neglected; for 
systems in which there is more than one component in the gas-phase, 1/Kg 


is generally small enough to be neglected when compared with H/K The 


1 


possible “exception is the case in which K. is very large because of a 


L 


fast reaction in the liquid phase. 


° = Me 4 (e) 
Values of K are often given as enhancement factor defined as KK, 


where ae is the physical mass transfer coefficient or in case of slow 
reaction regime Kk, = eae Thus the determination of K, generally depends 


on the evaluation of oe 
; ’ oa 
A variety of empirical correlations have been proposed for K ie 
bubbling systems, without any mechanical agitation. Calderbank and 


Moo-Young (51), Griffith (52), Johnson and Akheta (53), Miller (54), 


Soo (55) and Teller (56) have proposed the following type of correlation:- 


* 


saiteh oF Sdesigin et fosorgas scasteiess 
~esy to emiet ni pX .tnsioiitsos sstensiy Seen seedq-268 
9.1 ,csonsteles2 sesdq-binpil bas 2 
‘ Co Bal | 

(8 -&)  « + * @ * . @ 6) 8 © * # * * + * * z =a 4 = Pe 
| oe Se ie 
 ineisitises noitatoeds Isvimeto eesdq-biupit Isubivibat = orate” 
ah} 
emis\digaol fi 
: ; a > eal 4 
\eiom ,taetsitises sstensc7 cesm szesdq-2sg Leubivibnst = . 
Ss eeetg-smit-ssis ; 
: ao 
eyug B sotiw al meraye ooitose: bivpil-asg “sqyt-slddyd! yas al 
3O% gbetoalgen od aso .4\1 .sonstcieoy Sapig-esg siz , bsyolqms * 62 - 
\ 
gi\i .eendq-263 edt ai tascoqmes amo aedt sim at ersdt doinw of a 7 
ant INH Aatiw bevegnos oeiiw: Rady ALBPe ad oF dguvoms ifsme Waeh al 
.sasriq biupil et at coitonen jon} = 
“Pa? za param iotSet Tnomesasdas 25 nevis aatio om 1 te cau =) 


& To seusosd syisl yiev =2f 7 doinw uk sess sit wi NOLTGsoxs & 


wole to ease se fo Sasioitieos tiencit easm serio sit ak ‘2 a f 

7 

ebasgst) ti ehey: gi 70 noite tpikarteteb od> 2uhT op = 2? “To —_ 
vie 

| Pa 6 ee 


Fs svsii eared tind >ixkgms 


— 


46 


Oo 
K Vv 
oe ee Gis 
= av) rae e ° © © © e e e Cy © © ° 6 © . e e 6 
UL, it L D 


where U,, and L are characteristic of the liquid velocity and a 


dimension, respectively; a and a, are constants. A literature 


ee: 3 
review on various equations is given by Valentin (31). 

A correlation proposed by Hughmark (57) appears to be useful. This 
correlation applies to single bubbles an liquids, for liquid drops: in 
liquids, and for single spherical surfaces exposed to flowing air or 
liquid streams. The predictions are claimed to deviate by 15% from 


experimental data. The general form of the correlation for liquid-phase 


mass-transfer coefficients is presented as: 


b 
0.484 0.339 ,“D® 0.072) k <6. 9) 


where, D = diffusivity 


qd. = bubble diameter 


ays be are constants 
for single bubbles aint Oi OGLs, by = 1.61 
for bubble swarms ay = Oe Ode (ee send 


k 
For bubble swarms the velocity in the Reynolds number is represented 
by the slip velocity between the bubbles and the liquid. 


For cases where 1/3 


are >>2, 


the liquid-phase mass-transfer coefficients for single bubbles (SB) 
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and bubble swarms (BS) are related as follows: 


(K,°)Bs (a, )BS 
Sree Ss Ws Sl 6 al sen eUe) Belews ee Bue CSO) 


(K,°)sB (a, )SB 
Hughmark's correlation thus indicates that mass transfer 
coefficients for single bubbles are greater than those for bubble 
Swarms. 
Lochiel and Calderbank (58) present theoretical equations for 
calculating mass transfer from single bubbles in gas-absorption. 


; : 5 fe) 
According to them, when the continuous phase is pure water at 20 C, 
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which is the solution of Boussinesgq equation (58) for transfer around 


spheres in potential flow. 
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Calderbank and Moo-Young (51) have suggested the following equation 


for bubble swarms of average bubble diameter < 0.25 cm, 
ELS pe he} 
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and for average bubble diameter >0.25 cm, 
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where , oF density of continuous phase, gm/cc 
Pa = density of dispersed phase, gm/cc 
u_ = viscosity of continuous phase, gm/cm. sec. 


g = acceleration due to gravity, on/seeee 
D = diffusivity , om” /sec. 


Sharma et al. (50) have given a graphical correlation for Ka vs 


Use using various two-phase systems. Braulick et al. (60) have given 


some interpretation of literature data enabling one to calculate Ka 


corresponding to any superficial gas velocity. 


Most of the correlations for K have been developed for the air- 


water system. Extension of these equations to other systems is not 
definite. The conversion from one system to another could be made by 
the following relationship developed by the AIChE Distillation 


Committee (60). 
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Ce Gas Hold-up :- 


The gas hold-up in a sparged contactor is important in 
determining residence time and interfacial area for mass transfer 
and has generally been determined by directly measuring the height of 
aerated liquid and that of clear liquid without aeration. Thus the 


average fractional gas hold-up is given as 
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fractional gas hold-up 
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£ the expanded (aerated) liquid height 
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Yoshida and Akita (61) have Tene experimentally the effects 
of various parameters on gas hold-up, a summary of which is given below. 
Effect of Electrolytes on Gas Hold-up:- 

Yoshida and Akita (61) found that the fractional gas hold-up in 
water and electrolyte solutions in a 7.7 cm diameter column was slightly 
larger than in columns having a larger diameter, probably due to wall 
effects. This observation is in agreement with that of Fair et al. 
(62). Yoshida and Akita also observed that gas bubbles in water were 
generally less than 10 mm in size, whereas bubbles in electrolyte 
solutions consisted mostly of very fine bubbles and were dotted with 
larger ones of several millimeter diameter. The occurrence of small 
bubbles in electrolyte solutions can be explained by the 
electrostatic potential at the gas-liquid interface. 

Biteat or Gas. Rate :- 

Gas hold-up in sparged contactors varies directly with superficial 
gas velocity for the air-water system and for the air-sodium sulfite 
system. Graphs of gas hold-up vs superficial gas velocity have been 
given by Yoshida and Akita (61), Fair et al. (62) and Hughmark (57). 
Effect of Nozzle Diameter:- 

According to Yoshida and Akita (61), fractional gas hold-up is 
independent of nozzle diameter (range: 0.2 to 0.4 cm.). This can be 


explained by the fact that, in the range of gas rates studied in 
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their work, gas flows out of the nozzle as a continuous jet stream 
and is then split into bubbles by the turbulent motion of liquid in 
a zone several inches above the nozzle. 

Effect of Liquid Height:- 

Yoshida and Akita (61) found that fractional gas hold-up did not 
depend on the liquid height when the column height was over one 
meter. Since conditions in the vicinity of the gas inlet nozzle are 
different from other parts of the column, some end-effects are 
expected in shorter columns. This observation agrees with results 
given by Fair (12), according to whom the fractional gas hold-up 
increases as length-diameter ratio (2, /4,) is decreased, but is 
practically independent for Z,/d.. at 10 and above 10. 

Effect of Temperature :- 

Fractional gas hold-up is unaffected by temperature in the range 
ar | storage 0 (61). 

Hold-up Determination from Empirical Correlations :- 

Koide et al. (44) have suggested the following correlation for 
estimating fractional gas hold-up. The slip velocity, U. for 
bubble swarms is defined by 
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where Ur, a linear liquid or gas velocity. The ratio of bubble slip 
° 
velocity to terminal velocity, ie correlated with the liquid 


properties for two ranges of average bubble diameter: 
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PE) 
Ue Views 05,07 O78 CL ee) 


\ oto | was 
ai blupil %6 Aoisom tasludws ods yd selddudyoml tilge sds ab | 5 
. dt eveds - fgat Sya 
aissonm 9 eee +o re 


msarte taf aounitnod s 25 sissod sit Yo wd 


— 
> 


-:F gtsh b pid 2. oe : 
ton bib qu-biod esy, feacitosit tsd7 bavol (£8) s2idA bas sbifleoY 


‘ ' ar 
ano weve 26w tdgied.amiloo ers venwW tigishBiuplt ent 1 ie ; s" 
stata 7 


sxe efostie-bas smoe ,ammioo adi 26 etheq Tere mort 1RTA « 


; Fy 7 a . _—*- rio F ss 
eve efSson telat 2eg siz to vrinisiv edz af enengibage aoni2 


' 
® o . 
* 


aslwebs fatw 2e6e%es noltsvneedo <1nT -comifos yetsore mf barseqne 

. ae : 

* tw-biod: as Sencitostt sd° moun of gaibreoss <(Si) tied vd asvig i 

si ted’, bedsstosb <i (56N 8) oltser tetomp rh-dtahel es acepetont . yr 
; iy , 

of svods bn Ci te ,b\,8 162 Masbnsgebar gtaorzammg | 5 4) 


. n@T 


; 


. 
-1BautEIS 


sgus1 Sd? al stuteisqmes vi bsissiisay et qu-blod eeg lenoizoset , Fr 
hay. 9%0e es SE He: a 
-:enotisisa10) Iscinigrd mort fens. £2 rtd, , to - 7 
tot noltsletso> griwolloi adit iat-ackanan eved (#4) .is ts abion | 
| et 


Nol .U .ytioofsy gife sil .qu-biod gey Isqoitosd? parsers 


a > 


yd bemites ai anrrawe atddird 7 
U U 
Ot is td . a) (48 3 oe * ye * . a * -_ * . _ © * & a — awe fi 
(vt.@) . Ss Sew ¥ 
ae 


(ile efddvd 46 oltea Sat «ytinclev 2esy a0 -biup tt ‘eentl. * 5, U onedw 
Biuprt sdt aalw bassisyv102 esw me | CL ee  Sietigte yris 


7 Cad . 


= 


a 7 
on he Wea 
: as ; . i 3 yaa * 
tute a 


i 
> 


s 7 = 7 7 7 
ae TS ne es eee 


e 
_ 


2 
qd gp 


= Be me OP ANOD See ete tee a ie ge ule aw gee ee a) 


Large Bubbles:- d > 0.4 cm. 
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The terminal velocity may be determined experimentally, obtained from 
the literature (63) or estimated as shown by Perry (41). 

The average fractional gas hold-up of bubble swarms can now be 
estimated if negligible bubble coalescence is assumed. The procedure 
LexeteGa ie Eind qd. using an appropriate correlation, (b) determine a 
value of Ma based on ds. and (c) estimate ¢« by trial and error from 
Squationss(S.7)i and (5.18) or (5.17) and (5.19). Koide et al. claim 
an accuracy of +30%. 

Another recent correlation and probably the most useful is that 
of Hughmark (57). He obtained data in the bubble regime for air 
with water, a sodium sulfate solution, kerosine and a light oil ina 
l-inch tube. Gas hold-up was found to be correlated as a function of 
the superficial gas velocity for the air-water system at zero liquid 
flew. «The correlation applies to cocurrent liquid systems if the hold- 


up is defined by 
Un 
pty aotaguaphing. ta Adperaions: wet KoRBy, 2/5/20) 
S) 
The Hughmark data were used to evaluate the effect of liquid 


physical properties. The data indicate that hold-up for these systems 


can be correlated with the term 
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of liquid, which reduces to U 


pes) 
) / » where o = surface tension and P= density 


Usa 


for air-water system. 


SG 
From Highmark's correlation for e¢' vs Ugg (a EOS, e can be 
L 
back-calculated by trial and error for any particular Use anc Users 
using the following relation (57): 
U U 
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Gas hold-up data of other workers e.g. Fair et al. (62), Yoshida 
and Akita (61) and Towell et al. (64) compare pretty well with this 
correlation within a maximum absolute deviation of 21.5%. All these 
data were taken at 1 atm. pressure. However, the Neusen data (65) for 
steam-water at 600 psi in 2.9 inch pipe show an average absolute 
deviation of 32.5%. 
D. Interfacial Area: 

The first attempt to measure interfacial area was made only in 

1955 (66) and to date some independent techniques have been developed. 
Measurement techniques may be broadly divided into three different 
categories namely (a) optical methods, (b) photographic methods and 
(c) analytical-cum experimental methods. The details of optical methods 
can be found in the reference (67). Only photographic techniques and 
empirical correlations will be discussed here. 
Photographic Techniques :- 

This involves actually photographing the dispersions and then 


measuring the size of disperison by either optical methods (e.g. a 


a2 


cathetometer or a microscope) or other precision instruments. In addition, 


the number of spherical particles per unit volume of the mixture may 


‘tide = = 8 Sas a sostaue = 
ed nso 3 Nee, U ey "4 yor aotted iain 


CES, 20 iw! ee et okt (- a 
SbisiaoY , (8) ip. te tist .g.8 exexcow sades ites 
tds diiw Lew yttexq exeqnon (48) .Is to Lhawe? baw (£8) 
neat IIA «#2.18 to Aolwsiveb siwloeds ——e 5 aidtiw noltelers - wi 
in (28) eteb mediel odt ,»devewcl! . .ssvetesag ssueeoiq .mae, 1 30 cols sew stab 
stiloeds eybzsve. n= worle sqiq doad OLS AL reg O08 Fale ‘ “aa te ws 
ie &t te wyatt se i 
ab yino sham asw bets Isicsiretal aiwessm ot daniel terit sdT ) 
-bacolaysb weed syed asupiniss? tnebcaqsbat suse eteb et bas (28) 
Jgekeitibe ssant oral bebivib yibseoud od vem 2stipintoes t ‘ 
Bas: ehodsam -oidesrgesedg (d) ,abcodtem ine tage (ei) Yemen est 
ebedtem Issistqo, to elisteb si? .aborisem Levaeniaeque mo-Laoity . 
baa esupiatises sligeigeroiig ylad .( tS) sonsistex sit af 
Sisi bseavocib a eteaasaene 
san» cea 


‘ads oy eaoberedeie adt snide CERy:| 


23 


be actually counted, once the hold-up of gas is determined by any 
of the previously described methods. A mean bubble diameter called 


the volume surface mean diameter may be calculated by the following 


equation: 3 
yi 
pieay cnt 
Ve ~ 2 
aecell 
ae 
and then a! = be/di. gives the interfacial area per unit volume of 


mixture where € is the fractional volume gas hold-up. 

Sharma et al. (59) recently studied the effects of various 
parameters om interfacial area in a variety of gas-liquid systems in 
bubble-columns. The values of 'a' were evaluated by using the 
theory of absorption accompanied by fast pseudo-first order reaction. 
A summary of the effects observed by them is presented below. 


Effect of Z,/d,. Ratio:- 


For a column diameter of 6.6 om and superficial gas velocity = 
21 cm/sec, effective interfacial area was unaffected when Z,/4, was 
varied from 4 to 12. Similar observations have been reported by other 
workers (61,62,64). 
Effect of: Column Diameter:- 

The values of effective interfacial area obtained at the same 
superficial velocity of gas were found to be practically the same in 
the range of their (59) experimental conditions which were: superficial 
gas velocity = 15 om /sec to 40 cm/sec and column diameters = 6.6 cm 
te.3se7S%em. 
Effect of the Type of Distributor:- 

The values of effective interfacial area obtained by using different 


designs of spargers (59) (single-tube distributor, cross-type 
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distributor, ring-type distributor) in a 9 cm- diameter column were 
found to be the same when system properties and hydrodynamic 
conditions were unchanged. Similar observations were reported by 
other workers (68). 

Effect of Gas Velocity and Gas Density:- 

Sharma et al. (59) found that the effective interfacial area 
varies as the 0.7 power of the superficial gas velocity. 

When CO, diluted with Hy» air, or Freon-12 was absorbed in 
aqueous NaOH solution in the 6.6 om. diameter column, at constant 
superficial gas velocity, the effective interfacial area, a, was 
found to remain physically constant. This indicates that the change 
in the density of the gas is unlikely to affect 'a'. 

Effect of Pressure :- 

Tadao and Takeya (69), who studied the hydro-dynamics in a gas 
bubble column under pressures of 1 atm. to 6 atm., found no effect 
of pressure on bubble size or bubble rise velocity. Thus at a given 
superficial velocity of the gas, when the column is operated under 
pressure, it is very probable that the values of 'a' will be 
practically the same as in a column working at atmospheric pressure 
for equal values of the superficial velocity of the gas. 

Effect of Ionic Strength and Nature of Ions:- 

Yoshida and Akita (61), Braulick et al. (60) and Marruci and 
Nicodemo (70) have stressed the importance of the effects of ionic 
strength and nature of ions on gas hold-up and gas-liquid 


interfacial area. Sharma et al. (59) show that with an increase 
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in both viscosity and the ionic strength of the solution, the 
effective interfacial area increases considerably. When the 
viscosities of the solution are the same, aqueous non-electrolyte 
solutions give a value of effective interfacial area much lower 
than that in the case of electrolyte solutions. 

Effect of Non-aqueous Solvents :- 

It was found (59) that at equal viscosities the effective 
interfacial area increased from 1.7 em? CORA ewl se om + when the 
solvent was water and isopropanol respectively (surface tension: 
22 dyne/com). A decrease in surface tension thus leads to an 
increase in the effective interfacial area. 

Empirical Correlations :- 

Using an extensive collection of experimental data from tanks 
in which gas was dispersed by an agitator, Calderbank (71) obtained 
the interfacial area for mass transfer:- 
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ate 
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where Po agitator shaft power during dispersion 


ie) 
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as volume of agitated mass 
ye = superficial gas velocity. 
Vi. = terminal velocity of bubble 


a = interfacial area per unit volume of liquid-phase. 
The shaft power per unit volume of agitated mass can be estimated 


by the following equation given by Towell et al. (64). 
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CHAP TBR: Vil 


EXPERIMENTAL 


A. Kinetic Measurements :- 

This section describes the experimental equipment and operating 
procedures used to study the kinetics of the oxidation of acetaldehyde 
to acetic acid. 

(a) Experimental Equipment :- Figure 1 is a schematic of the equipment 
used. 

The reactor used in this study was a Parr-high-pressure reactor 
which had a capacity of 2000 ml. It was constructed of carpenter 20 
stainless steel and the pressure and temperature limits were 2000 psig 
and 350°C respectively. The reactor was equipped with an adjustable 
speed stirrer and stirring speeds up to 3500 rpm could be attained. An 
electric heater was used to heat the bomb. Water was circulated through 
a coil in the reactor to provide cooling. Temperature was measured by 
an iron-constantan thermocouple (type J) and controlled by (1) an on-off 
automatic temperature controller actuated by the controlling thermocouple 
and (2) by manually adjusting the cooling water flow rate. A 
refrigeration unit was used to produce cooling water (temperature ihe 
for the peactor. 

Charcoal and calcium chloride absorbers - 2"0.D. x 18" high lucite 
cylinders - were used to remove moisture and oily substances from 
oxygen to be used for the reaction. Glass wool was used to prevent 
entrainment of particles. 


The reactor pressure was maintained constant by means of a 
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pressure-regulator in the oxygen-inlet line. A one-way valve prevented 
the back-flow of liquid reactants from the reactor into the gas inlet 
line. Heat from packing glands was removed by cooling water circulating 
around them. 
Sampling of the reaction mixture could be accomplished by means 
of a sampling tube built in the reactor. Any leaks through the packing 
glands (in the stirrer assembly) could be detected by means of a gas- 
leak detector plastic tube immersed in water. During or after the 
reaction, as the need may be, the gases in the reactor could be purged 
through the gas-outlet valve. The gas-outlet line was connected through 
tygon tubing to the exhaust hood to prevent dispersion of chemicals in 
the laboratory. 
(b) Experimental Procedure:- The following steps were followed to carry 
out the reaction in the stirred-tank reactor. 
1. The reactor was cleaned thoroughly by using distilled water 
and dried with air. 
2. The system was checked for leaks. 
38. The pressure regulator, stirrer-speed controller and the 
temperature controller were set at the pre-determined values. 
4, The heater was brought to temperature. 
5. The bomb was charged with measured volumes of acetaldehyde 
and acetic acid containing a known weight of dissolved 
catalyst. 
6. The bomb was closed quickly and placed in the heater. 
7. Oxygen gas-flow was started, purging out the air from the top 
of the reactor and pressure in the reactor was raised to a 


fixed value, keeping in mind that any further temperature 
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increase would increase the pressure; the gas outlet valve 
was then closed. 

As soon as or just a little before the temperature and 
pressure in the reactor had been achieved, cooling water 
lines were opened and the stirrer was started. 

The pressure gauge on the reactor had to be watched 
continuously and the pressure readjusted by releasing some 
gas through the gas-outlet. 

Liquid samples were taken initially and at intervals of 5, 
10, or 20 minutes as the situation warranted. 

The sampling test-tubes were stored in the deep-freeze so 
that they’ were° cold prior to use. 

Before collecting any sample in the test-tube, a little bit 
of the first part of the liquid coming out of the sampling 
tube was discarded so that the sample collected was 
representative of the liquid inside the reactor at that 
time. 

Since the boiling point of acetaldehyde is eal eh the samples 
in test-tubes were quenched immediately by placing them in 
an ice-bath or deep-freeze. The temperature of the liquid 
samples’ was reduced to about o°c within 2-3 minutes. 

The reaction was run until the conversion of acetaldehyde 
was almost complete, which had to be determined by analyzing 
the liquid samples, consequently sample analysis was carried 


out while the run was in progress. 
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15. At the conclusion of the reaction the heater and the stirrer 
were shut down, the bomb was cooled, the gases in the bomb 


were released to the exhaust and the bomb was cleaned. 


(c) Experimental Conditions Used in the Stirred-Tank Reactor:- 


On the basis of the literature review and the features of 
experimental apparatus, the following range of experimental conditions 


were used. 


Pressure, psia ARS es Shar, KY, 7 
Temperature , ae OP HO, 50 
Catalyst (Mn (OAc), 80 > 250 
concentration, ppm. 

Volume % of AcH 10 d 20 
0D aH? 

Ree tls 500 Za 3000 
Volume of Liquid 500 ; 1000 


Mexture | ces. 


Manganese(ous) acetate (powdered) was used as a catalyst for this 
reaction because it gives the highest efficiency for the production of 


acetic acid, 


w 


Sparged Reactor Study 


In this section the equipment and the operating procedures used 
for the sparged reactor study are described. 
(a) Experimental Equipment :- Figure 2 is a schematic of the bubble 
reactor and associated equipment. 


The bubble column was 4 ft. high and was constructed of four-inch 
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I.D., Q.V.F. pyrex glass (72) sections and was capable of withstanding 
up to 40 psig working pressure. The bottom section 30 of the column 

was a 6" long x 4" I1.D. 316-stainless steel section. An orifice plate 
of any desired design could be inserted at the bottom of the column. 

The second section 35 of the colum was also a 4" I,D. x 6" long 
stainless steel 316 pipe, which had two pyrex glass windows through 
which the bubble streams at the orifice tips could be observed and 
photographed. The temperature of the column was controlled by recir- 
culating the reaction mixture at high rate through a heat exchanger 
where heat could be added cor removed as required. Liquid was withdrawn 
from the top of the column at the point 14 through a 3/4" flexible 316 
stainless steel pipe and through a gate valve to a surge tank 15 where 
the liquid and any entrained gas-bubbles were separated. This vapor 
returned to the column through a line. The recycle line from the surge 
tank was connected via a 3/4" flexible 316 stainless steel pipe to a 
centrifugal pump 16. Flexible stainless steel pipes were used to 
minimize the transmission of vibrations of the pump to the glass column. 
The recycle liquid was pumped through a coil and shell-type heat 
exchanger 18, where it was heated or cooled depending on the requirements. 
It then flowed through a rotameter and an one-way valve to the liquid 
entrance line 12 near the bottom of the column. 

Three thermocouples (Iron-constantan, type J) 32 were inserted in 
the column’ through 1/2" thick stainless steel 316 rings placed between 
each section of the column. Cold water was used for cooling the recycle 
liquid whereas hot water was used for heating. The hot water was 


obtained by heating cold water with steam in a shell and tube-type heat 
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exchanger 19. A solenoid valve 21 was installed in the steam line and 
was connected to an automatic on-off temperature controller 36 actuated 
by a temperature probe 33 inserted in the column. The solenoid valve 
opened when the temperature in the column was lower than the set value. 
Thus a constant temperature could be maintained in the column by heating 
or cooling the recycle liquid. The flow rate of water and the steam 
pressure had to be adjusted manually. The recycle liquid flow rate was 
about 2 gals/min. 

Oxygen gas at a constant pressure flowed from a gas cylinder 
through the needle valve 1 and similarly nitrogen gas flowed through 
the needle valve 2. Both gases were passed through a pressure-regulator 
to an activated charcoal absorber and to an anhydrous calcium chloride- 
absorber, where oily substances and moisture from the gas were removed, 
respectively. Dry and oil-free gas then passed either through rotameter 
3 or 4 and through a pressure gauge, through a needle valve 5 and 
through’ a one-way valve to the orifice plate 13 in the bottom of the 
column. The one-way valve prevented any back-flow of liquids from the 
reactor into the gas-inlet line. All gas-inlet lines from needle 
valves 1 to 5 were 1/4" copper lines and from needle valve 5 to the 
column entrance, they were 1/4" 316 stainless steel. 

The equipment used was capable of handling continuous flow of 
both gas and liquid. However, in this study the system used was 
"continuous gas but batch liquid'. Glacial acetic acid (99.99%) was 
stored in a polyethene carboy 9 and was charged to the column with a 
centrifugal pump: A single orifice’ 29, 0.095" 0.D. (0.078" I.D.) x 


4"" long constructed of 316 stainless steel tube was used. The gas- 
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phase flowed in the column through the orifice at 13 and liquid-phase 
flowed through a side entrance at 12 and through a number of holes in 
the top surface of the section 30. 

A pressure gauge (0 - 100 psi) 27 was installed in the top section 
of the column. The pressure gauge could be screwed out from the tube 
31 and acetaldehyde was introduced into the reactor through this opening. 
A pressure-relief valve 26 was installed at the top of the column and 
could be set at 30-40 psig. A shell and coil type condenser 25 (stainless 
steel 316) was employed to condense vapors in the exit gases. Cold water 
was used for condensing the vapors. A drain valve 6 was installed at 
the bottom of the column. Liquid samples were withdrawn through a needle 
Weave AJ 3 
(b) Measurement of Gas Hold-up:- 

Gas hold-up measurements were made for the oxygen-acetic acid 
system at 20 psig and Ho°e. Gpneee were the conditions at which the 
reaction was carried) out in thé column. 

Operating Procedure: - 
The column was filled with acetic acid to a height of about 40 
inches. Oxygen was obtained by opening the valve 1 and the flow rate 


was recorded with one of the rotameters 3 or 4. Liquid levels at 


ip) 


various gas flow rates were measured using a cathetometer. Measurement 
were made with the recycle pump in operation. 
(c) Measurement of Bubble Diameters: - 

The bubble diameters for various gas flow rates were measured 


by a direct photographic technique. A high speed (1/1000 sec) camera 
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was used with a flood light at the back of the column. A concavO-plano 
lens was designed of lucite with water in the lens to fit one side of 

the column, through which preliminary photographs of bubbles were taken. 
The lens was expected to reduce the distortion of bubbles in the 
photographs, but significant distortion did not occur in the absence 

of the lens. The inside diameter of the column was determined by placing 
a scale horizontally in the centre of a section of the column with water 
inside it and taking photographs from various distances. From four 
different photographs taken from varying distances, the following 


relation was calculated. 


Diameter of the column in the photograph _ Length of one inch scale division 
Actual inside diameter of the column in the photograph 


one inch 
The length of both the l-inch scale division and the diameter of 
the column in the photographs were measured and the actual inside diameter 


of the column (x) could be calculated from the above relation, i.e. 


a length of l-inch in the photograph 


diameter of the column in the photograph in inches 
From the four different photographs mentioned above, it was found 
thay x =o 2 evant oles, WO" Ssrs0"Gand hence an average yalue of x = 
3.45" was used=as the actual inside diameter of the column. This also 
proved that there was not any appreciable distortion when photographs 
were taken from varying distances. The following relation was used 
for calculating the actual diameter of each bubble. 


x(bubble diameter in the photograph) 


Actual diameter of a bubble = diameter of the column in the photograph 


The bubble sizes were measured by means of a microscope. 
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(d) Measurement of Rate of Conversion of Acetaldehyde :- 


The reaction: oxidation of acetaldehyde to acetic acid was 
carried out in the sparged reactor at the following conditions: 20 pSsie, 
uo°c, 140 ppm catalyst concentration, constant liquid height but at 
two different oxygen flow rates. The initial concentration of 
acetaldehyde in the mixture of acetaldehyde and acetic acid was approx- 
imately 10% by volume and pure oxygen was used. 

The column was initially cleaned by charging it with acetic acid 
and passing nitrogen through it. The liquid was circulated so that 
any water that might have remained in the column was dissolved in the 
acetic acid. The column was charged with acetic acid and drained 2-3 
times to make sure that there was no water in the system. 

The following steps were taken to carry out the actual reaction 
run :- 

1. The system was charged with acetic acid. The condenser 
cooling water inlet valve 24 was opened. A predetermined 
volume of pure acetaldehyde (99.98%+), which was stored at 
about O°F in a deep-freeze, was then poured in the column. 
The recycle pump was started. 

2. Nitrogen gas-flow was then started and its flow-rate was 
adjusted to a known value. Pressure in the column was raised 
toe20ipsig, 

3. The temperature-control system was switched on, and the 


P , fo) 
temperature in the column was raised to WORE. 
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4, Nitrogen-flow was stopped and oxygen flow was started. Since 
the reaction is exothermic, the temperature in the column has 
to be re-adjusted to 4O-C by adjusting the cooling water valve 
23 and the steam-pressure-regulator. All these Steps on 
bringing the system back to 40°C and 20 psig had to be done 
as quickly as possible. 

S. Samples of liquid were withdrawn initially and after intervals 
of 5, 10 or 20 minutes, as required and were quenched in a 
deep-freeze, 

6, Analysis of liquid samples were done as soon as possible. 

After the completion of the reaction, oxygen flow was stopped and the 
temperature-control system was switched off. The system was drained 


after the liquid in the column was back at room temperature. 


cs Analysis of Liquid Components:- 


The liquid mixture contained mainly acetaldehyde and acetic 
acid and had to be analyzed for acetic acid content. This was done 
simply by titrating the liquid sample against standard sodium hydroxide 
solution: 

CH ,COOH + NaOH > CH,COONa sr H,,0 

Phenolphthalein was used as an indicator. A magnetic stirrer was used 
to stir the liquid. Titration was done in an inert atmosphere of 
nitrogen to reduce the possibility of any further oxidation of acetaldehyde 
while titrating. However, if the titration was done quickly enough, it 


was found that there wasn't any appreciable difference in the results 


when titration was done in an inert atmosphere or an ordinary atmosphere. 
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It may be noted that in the presence of dilute sodium hydroxide 
solution (10%) which acts only as a catalyst, aldol condensation would 
take place between the aldehyde molecules. However, it will not 
interfere with the titration because sodium hydroxide does not react 


with acetaldehyde but only acts as a catalyst. The aldol condensation 


reaction is shown below: 


OH 
CHG. SS 0a =Cc 
ae O + HCH,CHO % Giger 
OH Aidol 
-H,0 
VW Z 


CH CH = CHCHO 
(Crotonaldehyde ) 


This was checked by titrating 1 ce of acetic acid with N/1 sodium hydroxide 
solution and then 1 cc of acetic acid + 1 ec of acetaldehyde mixture 
against N/1 sodium hydrcxide, In both cases, sodium hydroxide consumed was 
the same. Thus ecetic acid was analyzed by titration in this work. 
phalysis for Acetic Anhydride, Peracetic Acid, Water andor olher by-prodicis: 
As was seen in section 2, depending on the experimental conditions, 
ene. suemperature, Catalyst, purity Of reactants, theres could be) products 
G@iher tian acetic acid in the oxidation of “acetaldehyde to acetic acid: 
To make sure that the products did not contain any appreciable amounts 
of these by-products, occasional checks were made on the samples by gas- 
chromatograph and mass spectrometer. In the products formed in various 
reaction runs, no traces of peracetic acid or acetic anhydride were 
Pood In some cases, only traces of water were found but the percentage 
of water was less than 0.1% and thus within the limits of experimental 
error, it was assumed that there was complete conversion of acetaldehyde 


to acetic acid. 
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Gas-Chromat ograph :- 

Poropack-S column, 6' x 1/8", was employed at 120-180°¢ using 
helium (carrier gas) at the rate of 30 c.cs/min for analyzing the 
reaction mixture in the gas-chromatograph. The elution times at 150°C 
for water, acetaldehyde and acetic acid peaks were 0.4 min, 0.6 min 
and 3.2 minutes, respectively. Mass spectrometer was used to detect 
any acetic anhydride content. 

Chemicals: 

Acetaldehyde and acetic acid were supplied by Chemcell Limited, 
Edmonton and were reported to be over 99.9% pure. The catalyst, 
manganese(ous) acetate, powdered, was purchased from Fisher Scientific 


Company. 
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CHAPTERe VIL 
RESULTS 
Pe Kinetic Studies:- 
Kinetic data for the oxidation of acetaldehyde to acetic acid 
were measured in the stirred-tank gas-liquid contactor described in 
section VIA. The effect of the following variables on the reaction rate 


were investigated: - 


x Concentration of acetaldehyde 
Zs Concentration of oxygen 

On Temperature 

4. Catalyst concentration. 


(a) Determination of Operating Conditions for the Kinetic Regime:- 


Measurements were made at different stirring speeds to 
determine operating conditions in which the kinetic regime theory could 
be applied. In this regime the liquid must be physically saturated with 
oxygen and thus the reaction rate should be independent of the degree 
Of agitation. In addition the rate of reaction per unit volume of liquid 
should be independent of the volume of liquid used. 

Stirring speeds ranging from 600 to 3000 rpm were used. The 
reaction conditions for initial runs were: 

Catalyst concentration = 140 ppm 

Initial acetaldehyde concentration = 10 vol % approximately 
Initial volume of liquid = 500 ccs. 

Pressure = 24.7 psia. 


Data are shown in appendix B, table B.1 and the results are shown 


on figure VII.1. It is seen from this figure that the conversion of 
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acetaldehyde to acetic acid increases as the stirring speed is increased 
from 600 to 1600 rpm but then decreases at higher stirring speeds. 

The decreasing reaction rate observed at high stirring speeds can 
be explained by assuming that vortex-formation was taking place causing 
a decrease in interfacial area. 

Additional runs were made at 34.7 psia and 44.7 psia. Data are 
shown in tables B.2 and B.3 and the results are shown on figures VII.2 
and VII.3, respectively. The reaction rates for the 24.7 psia and 
34.7 psia were highest at 1600 rpm but for 44.7 psia, it was highest at 
2135 rpm. It appears that the stirring effect on the liquid particles 
decreases at high pressures and a high stirring rate is required to 
achieve saturation of the liquid with oxygen. It was also observed that 
the reaction rate did not decrease appreciably as the stirring rate was 
increased from 2135 rpm to 2900 rpm. 

The highest concentration vs time curve was taken to represent 
the chemical reaction controlled regime. To further check this 
assumption runs were made using a larger volume of reaction mixture 
(i.e. 1000 ccs). One run was made at 1600 rpm and the other at 2135 
rpm, both at a pressure of 24.7 psia. Data are shown in table B.4 and 
the results are shown on figure VII.4. The reaction rate at 1600 rpm 
was the same for both the 500 ccs. and the 1000 ccs runs, thus supporting 
the conclusion that the reaction rate was independent of the degree of 
stirring at this condition. At 2135 rpm, this was not the case. Thus 
the reaction rate was independent of liquid volume at 1600 rpm and 24.7 
psia and the kinetic regime theory could be applied at the condition 


where the highest reaction rate was obtained. 
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Ube barect or Acetaldehyde Concentration:- 

The effect of acetaldehyde concentration on the reaction rate 
could be obtained from any of the acetic acid concentration vs time curves 
obtained by operating in the kinetic regime. Either a differential or an 
integral method of analysis could be used. 

(c) Effect of Oxygen Concentration:- 

The equilibrium concentration of oxygen-gas in the liquid 
mixture varies with pressure and therefore by carrying out the reaction 
at different pressures, the effect of oxygen concentration on the reaction 
rate could be studied. The results obtained at three different pressures: 
24.7 psia, 34.7 psia and 44.7 psia are shown by the highest curves on 


figures VII.1, VII.2 and VII.3 respectively. 


Vd) btrect ot Temperature :~ 


Since the chemical reaction 
the stirring rate required to achieve 
temperature should also be sufficient 
rates used at 50°C were also used for 

fo) fo) : 
eo cer Cland 40°C, bothwat 324.7 psia. 
results are shown on figure VII.5 for 


(e) 


rate is higher at higher temperature, 
the kinetic regime at a higher 

at a lower temperature. Thus stirring 
lower temperatures. Runs were made 


Data are given in table B.5 and 


30°C, 40°C and also 50°C. 


Eiake\Ciun Ors Catalyst Concentration:- 


Reaction runs were carried out at two other catalyst 


: fo) 
concentrations 83 ppm and 249 ppm, both at 34.7 psia and 40°C. 


Data are 


given in table B.6 and results are shown on figure VII.6. 


(£) Check of Reproducibility:- 


: fo) 
A run was repeated at 140 ppm catalyst concentration and 40°C, 
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34.7 psia and 1600 rpm. Data are given in table B.7 and the results are 
shown on figure VII.7 which also shows the run taken earlier at the same 
conditions. The reproducibility is seen to be good. 

(g) Catalyst ACTIVI ty t= 

A reaction run was made at Ho°C, 34,7 psia, 140 ppm catalyst 
concentration, 1600 rpm and initial acetaldehyde concentration of 20% 
by volume. Data are shown in table B.7' and the results are plotted on 
Figure VII.7' as concentration of acetic acid vs time. Acetic acid 
formed in this run was used for the next run. Fresh acetaldehyde was 
used and the make-up catalyst (since there was some loss of catalyst due 
to sampling) was added to bring the catalyst concentration to 140 ppm. 
The run was taken at exactly the same conditions as the previous run. 
Data are shown in table B.7' and the results are plotted on figure VII.7'. 
It can be seen from this figure that the 'concentration vs time' curve 
for the second run is lower than that for the first run indicating that 
the catalyst is deactivated after one run. 

(h) Determination of Rate Equation:- 

The review of literature in section II indicates that an 
assumption that the overall rate of reaction is first order in both 
oxygen concentration and acetaldehyde concentration may be justified. 

A rate expression, based on this assumption, was used to derive the 
following relationship between acetaldehyde concentration and time. 
The derivation, given in appendix A, takes into account the fact that 
the partial pressure of acetaldehyde decreases (partial pressure of 


oxygen increases) as the reaction proceeds. Acetic acid was supposed 
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to be the only product 


= 
P'-1.015x10 . i 
q x a Vee Da oH 
ray, ta ee remap ae ) = K,tP 
Pistols PA 
x10 (v, V,)b AcH AcH 


where b oH = concentration of acetaldehyde (volume %) at any time t 


4° ee aa 
AcH = tnitia value of D oH 


<= 
" 


vapor pressure of AcH at reaction temperature, psia 


<j 
' 


= vapor pressure of AcOH at reaction temperature, psia 


pf = (P-V,)5 psia 


-2 e) 
Pl=1 015x110 ~(V_=V.)b lea 
Let I' = ( a OS yee sle 
-2 fe) b 
ie 0 
IPM oll, OLS) 210) Ou V, db AcH AcH 
gyi so aie 


The assumed rate expression will be consistent with the data if the 
plot of InI' vs t.P' gives a straight line. The slope of this line then 
gives the value of Ke: 

Tables B.8, B.9 and B.10 show the processed data for runs at 50°C, 
140 ppm catalyst concentration and at three different pressures 24.7, 34.7 
and 44,7 psia, respectively. Figure VII.8 shows the plot of inI' vs 
ee for runs at 50 C, 140 ppm catalyst concentration and at all three 
pressures 24.7 psia, 34.7 psia and 44,7 psia. It is seen that the plot 
is a straight line and thus the assumed rate expression is consistent 
with the data. The slope of the line (K,) was calculated to be 7. 6leLO / 


(psia)(min) by a least square fit of first order. Because of the high 
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reaction rate at 44.7 psia and 50°C, data were integrated starting from 
95% (by volume) acetic acid, 

Tables B.11 and B.12 show the processed data from runs at 34.7 psia, 
140 ppm catalyst concentration but at uo°C and 30°C, respectively. InlI' 
vs tP' plots for the data in tables B.10, B.11 and B.12 are shown in 
figure VII.9, which give three different straight lines for three 
different temperatures. Slopes of these lines thus give the reaction 
rate constants at the corresponding temperatures:- 
From figure VII.9, 


= oon y Ly Peas MN. 


iB A 30°C 
: -3 ry Vee 
a nee = 3.46x10 L/psia min. 
~3 i ; 
com 50°C = 7-64x10 L/ps iar min. 


Arrhenius' Equation Fit for Temperature Dependency of Keim 


Arrhenius’ equation is 


o —-E/RT 
e 


i 


where, K a frequency factor 


R 


E = activation energy 
R = gas-law constant 
T = temperature. 


The above equation can easily be reduced to: 


fe) 
-Ink, =O F/R Ink, 


Thus, if -Ink, vs 1/T is plotted, it should give a straight line. 


The slope of the line = E/R. 


A plot of log Kp vs 1/T is shown in figure VII.10 and from cag Obed 


3 
plot E/R = 6.9x10°/°K. 
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. - E = activation energy 


= 13. 8x10° cals/g-mole 
Ke is calculated to be 1.48x10", 1/psia min, and the Arrhenius’ equation 


3 
can be written as K, Pee sio e 18+ 8x10 /RT 


Dependence on Catalyst Concentration :- 


Tables B.13 and B.14 show the processed data for the reactions at 
34.7 psia, 40°C and at two different catalyst concentrations ,83 ppm and 
249 ppm. LnI' vs tP' plots for these data and also for 140 ppm catalyst 
concentration are shown in figure VII.11. Straight lines are obtained 
for each catalyst concentration, which agrees with the assumed rate 


expression. Slopes of these lines (i.e. Kp) are given below: 


source Tae P,psia Catalyst Concentration ,ppm Kpoi/psia min 
83 2,02x10 ° 
Pies vileli) 40. 34.7 140 avuexlon 
249 7.00x10°° 


Figure Vil,l2éshows a plot of K. vs catalyst concentration on 4 


R 


)log-log plot and the slope of the line is almost one. Thus the reaction 
rate can be assumed first order in catalyst concentration. Figure VII.13 


shows a plot of K, vs catalyst concentration which gives a straight line 


R 
of slope (K,') = Soa 1/(psia min ppm). Hence the general rate 


expression is 3 
-13. 8x10 
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BiG: VIil.13 DEPENDENCE OF THE REACTION RATE ON CATALYST 
CONCENTRATION. 


Comparison Oar Experimental and the Predicted Results:- 


Gonsicer the Yun number 15.) P = 3h. 7 pera, T= 40°C, Se he 9.5%, 
\c 
catalyst = 140 ppm 
BiG ealor 
Predicted K, = 9.08x10' RI (eat. cone ) 
; dK, - 00x10" = Gat COnChs). Dp 


' 


=5 
2.92x10 (140 ppm) 


3. 46x10” 1/psia min. 


i 


Table B.15 shows the comparison of the predicted and the experimental 
data for t vs DoH and figure VII.14 shows the calculated and the 


experimental t vs. DoH plots, from which it can be seen that the rate 
expression presented here is quite satisfactory. 
Visual Observation: 

The reaction mixture which was colourless at the start of the 
reaction turned brown after the reaction started. This observation 
is in conformity with that of other workers (5,6). The color change 
noted is associated with a change in the valence of the catalyst from 


+++ 
@) 


ney to Mee (or Cane OMe ). If there is no color change, there 


is no subsequent reaction. 


Die Sparged Reactor Study: 
(a) Hold-up Measurements :- 


The oxygen-acetic acid system was considered to be 
representative of the reaction mixture and hold-up measurements were 
made for this system at uo°C and 20 psic. (A liquid cireulation rate co, 
fe gallons per minute was used. Results are given in table B.16 in 
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appendix B and are plotted on figure VII.15 as € vs Use (—— . = ) 
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FIG. VII.15 HOLD-UP CORRELATION FOR 0,,- AcOH SYSTEM AND COMPARISON 


WITH HUGHMARK'S CORRELATION. 
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Comparison With Hughmark's Correlation:- 


Hughmark's data take into account the effect of liquid physical 


properties on hold-up. The data indicated that hold-up for different 


systems can be correlated with the term U._(-= ; Tea Hughmark's 
Oo 
L 
data, however, were taken with zero liquid flow. From Hughmark's 
. 2 
coprelation for ¢€' vs U ( LES oS € can be back-calculated from 


SG fe) 


p 
L 
equation(5.19),as shown in section V.C, by trial and error for any 


particular ye and gr: Thus the data in the present work could be 


94 


compared to Hughmark's data. Table B.17 in the appendix shows the results 


62.4 ‘ j 
for Use ia7iate ay" vs € (experimental) and « (predicted from Hughmark's 
L 


correlation) and figure VII.15 shows the comparison of the two. It can 
be concluded that the experimental results are in good agreement with 
Hughmark's correlation. 

(b) Bubble Diameter and Interfacial Area:- 

Bubble Diameter: - 

It was noted visually that the bubbles formed in the oxygen- 
acetic acid system were much smaller than those formed in the air-water 
system. This demonstrated that the correlations available in the 
literature for the air-water system may not be suitable for the oxygen- 
acetic acid system. 

Bubble diameters were measured from photographs taken of the 
oxygen-acetic acid system at various oxygen flow rates at 20 psig and 
40°C. The oxygen flow rate was in the range of 47 ccs/sec to 334 ccs/ 
sec. The flow regime was determined from Govier's chart (30) and 
the system was found to be operating in the bubble regime. The bubble 


sizes were determined by scaling maximum (1) and minimum (w) dimensions 
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of the individual bubbles on the photographs. In general’, Moettor 
the observed bubbles were oblate spheroids; the diameter of a sphere 
equal in volume to the oblate spheroid, given by the following 
equation, was taken to be the diameter of a bubble: 

d= (ay)! * 


The volume surface mean diameter was calculated by the following 


equation 


Processed data for the bubble diameters are given in tables B.18 to 
B.22 in appendix B. The results are shown in table VII.1 and plotted 
in figure VII.16 as the superficial gas velocity of oxygen vs disieas 
The mean bubble diameters calculated from the empirical correlation 
(5.4) are also shown in table VII.1 for the comparison purposes. It 
must be noted that the experimental bubble diameters are much smaller 
than those calculated from equation (5.4). This was expected since the 
bubble diameter in electrolytes are smaller than in the air-water system. 
Interfacial Area:- 

The interfacial area in the bubble column is given in figure VII.17 
as a function of the oxygen flow rate. The values of interfacial area 


were calculated using the following equation: 
_ bE 
ae 

vs 


The processed data and the results are shown in table VII.1. 
Calderbank's equation (5.22) was also used for calculating the 
interfacial area and the values are given in table VII.1. They are 
smaller than the experimental values, demonstrating the fact that the 


empirical equations may only be suitable for the air-water system. 
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BUBBLE DIAMETERS, HOLD-UP AND INTERFACIAL AREA AS A FUNCTION OF SUPERFICIAL GAS VELOCI1 
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(c) Conversion of Acetaldehyde to Acetic Acid:- 
(ib) Experimental Results:- 


Two runs, one at an oxygen flow rate of 85 cces/sec and the 
other at 334 ccs/sec were made for the oxidation of acetaldehyde to 
eee eacid «at OC: 20 psig, 140 ppm catalyst concentration, 10% 
(volume) initial concentration of acetaldehyde and 40 inches liquid 
level in the bubble reactor. The data are shown in tables B.23 and 
B.24 in appendix B and the results are plotted in figure VII.18 as 
concentration of acetaldehyde vs time for both oxygen flow rates. 

(2) Modelling Sparged Reactor: - 

In order to select an appropriate model for the sparged 
reactor it was necessary to determine the reaction regime consistent 
with the experimental operating conditions. To do so, it was essential 
to estimate some parameters, i.e. Henry's constant (H) for the 
oxygen-acetic acid system, the diffusivity of oxygen in acetic acid(D) 
.) 


and the physical mass transfer coefficient (K For these purposes, 


L 


the reaction mixture was assumed to have properties of pure acetic acid. 
Solubility of Oxygen in Acetic Acid:- 

There are no data available for the solubility of oxygen in 
acetic acid. However, it was calculated from Prausnitz and Shair's 
equation (73) using the limited data available for the solubility of 
nitrogen in acetic acid. Calculations are shown in appendix A and the 
Henry's constant (H) for oxygen-acetic acid system at 20 psig and 
40°C was estimated to be 1660 atmospheres /g-mole 0,,/g-mole AcOH. 


Diffusivity of Oxygen in Acetic Acid:- 


Because of the lack of the data available in the literature for 
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the diffusivity of oxygen in acetic acid, it was calculated from 
Kuloor's equation (74) as shown in appendix A. It was estimated to 
be Da5H110” a ee at HOC 
Physical Mass Transfer Coefficient, Koi 

Calderbank and Moo-Young's (51) and Hughmark's (57) equations 
were used for estimating ee for the oxygen-acetic acid system. The 
results are shown in table VII.2. 
Determination of the Reaction Regime y= 

The diffusion time and the reaction time were calculated, as 
shown in appendix A, for the complete range of experimental conditions. 


The results are given below: 
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Thus ty<<tps indicating that the bubble reactor was operating in 
the slow reaction regime. 
Mathematical Model:- 

The following assumptions were made in order to analyze the 
reactor: 
(1) Liquid phase is well mixed and operated as a batch 
(2) Gas is passed through the liquid phase in the form of bubbles 
(3) Pure oxygen gas (99%+) is used as the gas phase 
(4) Gas phase resistance is negligible 
(5) Axial dispersion is neglected 
(6) Reactor volume is constant. 

fo) 


In the slow reaction regime, K, = K) and the general equation (4.1) 


for this regime, given in section 4, reduces to the following equation: 
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K Oo a! ( t dc. 
— (C'-C = = 
ee ») 1/2K.b , Co Fiancee we oe) na ee 


and the kinetic equation at a particular catalyst concentration is 


given by: 
db 
a AcH 
KEP a oH@o =e dt. e © © ® « o e ° ° e © e ° e e e ° ® eCie8) 


The above two equations were solved simultaneously by the 
Runge-Kutta method of fourth order to predict the performance of the 


reactor. The following initial conditions were used: 


The values of a', K, and e used in the model are given in table 


R 
VII,2. The best-fit° values for cae were determined by fitting the model 
to the experimental data for the acetaldehyde concentration profile 

in the sparged reactor. The computer program is given in appendix A and 
the model results are given in tables B.25 and B.26 in appendix B. 
Figure VI1.18 shows both’ the’ model prediction’ and’ the experimental 

curve for acetaldehyde concentration vs° time. The modei seems to fit 
the experimental data real good. The best-fit values for a apne” 0.015 
em/sec and 0.0065 cm/sec for oxygen flow rates of 85 ccs/sec and 334 
ces/sec, respectively. It can be seen from table Vil.2° that’ the values 
cag a predicted by Hughmark's equation for bubble swarms (57) and 
Calderbank and Moo=Young's* equation (51) for bubble swarms of average 


bubble diameter’ less than’0.25 cm are’ in reasonably good agreement with 


those required for the best-fit. 
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TABLES Vit v2 


PHYSICAL MASS TRANSFER COEFFICIENTS 


For 0, = 80 ecs/sec, £ = 0.0265; 
gs ecs of Iatad 
= 4,92 Se 
Ke ae g-mole AcH.min 
fe) s i) 
The best fit values: Ky a' = 0.0107 1/sec; a' = 0.7123 Ke = 0.015 om/sec. 
cm 


For 0, = 334 ces/sec, € = 0.112 


The best fit values: K at = 0,267 l/sec: a’ = 4.12 em’; x” = 0.0065 cm/sec. 


L L 
om 
85 cc/sec 334 cc/sec 85 ec/sec 334 cac/sec 
fe) fe) Oar OF 5 
Source K, Ko Ky .a KS <a 
em/sec em/sec 1/sec l/sec 
Calderbank's equation 
for qd. 2 Oy As) une OMOTRES On OUaS 0.0082 0.0474 
Calderbank's equation 
for d. > O525 Gi 0.0433 0.0433 0.0308 (5 ALAS} 
Hughmark's equation 
(1) for Single Bubbles ORO2Z 65 O22 On O85 ONO 775 
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a 
sose\mo 8£0.0 = op PRELN.O F'e pose\f VOLO,D = 1 
“a : ; 


.pee\m> 2900.0 = a, ae) GI.# = 's yaee\l Vos.0 = "5 Vv seoulav FF 
3 . 


‘Ms 
ase\bo BEE Dea\oo 2b agz\oo HEE sse\oo 88 
é © © a. 
‘ . | ' ‘ My A A : 
sen\ t sae\ f 5ae\ma oge2\no. 


See thee ee NRE A ee 


sottsupe 2 tas 
HTHO,0 © 2800.0 er£0.0 ALLO0 1 ALLO | 
. qo itsupe 2! 

BNL. 0 8060.0 SEHO.0 ceuo.0 mo 38.0 < 


_-nokaBape Re su 
a we 
avyo.0 #8AL0.0 850.0 2380.0  eekddyl afgade « 2 (L 
i » 7 v 


S10 .0 5300.0 acV00,0 T800.0 


103 


“dO LIvied TaApaVas SHEL NI 


CSA LON TN ) 
09 OS 


SENTRA pejoIpedg O 


Hee *TeLUsuTasdxy Vv 


cg *‘Tejueutuedxg gq 


SAIIAO¥d NOILVYLNAONOO GINOIT 


On o€ 0g 


sp a Os 


OT 


eR is: 


HOV 
AWNTOA ) q 


0 
(e) 


( 


104 


CHAPTER VIII 


DISCUSSION AND CONCLUSIONS: 

The performance of a two-phase sparged reactor can be estimated 
if the following information is available. 

(1) Description of the chemical reactions occuring and kinetic 
rate data. 

(2) A model which adequately describes the mixing pattern in 
the reactor. 

(3) Interfacial area. 

(4) Gas hold-up. 

(5) Physical mass transfer coefficients. 

For this study of the catalytic oxidation of acetaldehyde to 
acetic acid, since the reaction mechanism is complex, a simplified 
approach was used in that only the overall reaction was studied. A 
kinetic rate expression to describe the overall reaction was obtained 
which is first order in acetaldehyde - in agreement with the results 
of other workers who have studied the oxidation of aldehydes (4, ll, 
20, 21, 22), first order in catalyst concentration, which is also in 
agreement with studies made by other workers (2, 11, 16, 23) and first 
order in catalyst concentration, which is supported by the works of 
Bawn and Williamson (4) and Allen (9). The catalyst concentration used 
was in the industrial range (i.e. 80-250 Mn** ppm, equivalent to 
0.0253 - 0.079 wt % of Mn’ in acetaldehyde). 


The rate expression obtained in this work is given below: 


db 


3 
AO == 3 Li = hoe Baul: Ose AR 
dt 2S ACH a teas “DacytPo: (Cat. concn. , ppm) 
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where, DH = volume % of acetaldehyde in the liquid phase 


and Po 


rae partial pressure of 0, in the reactor, psia. 


2 

Venugopal et. al. (11) concluded that the catalytic oxidation 
of acetaldehyde to acetic acid was independent of the catalyst 
concentration which appears to be highly doubtful. Further, it is 
not at all clear what units are used in their rate expression (11) 
and there are inconsistencies in their paper, which make it impossible 
to use their rate expression even for comparison purposes. 

Studies in the four-inch bubble column provide information on 
bubble size, gas hold-up and interfacial area for the AcH-AcOH-0,, 
system. 

The fractional gas hold-up results are in excellent agreement 
with those predicted from Hughmark's graphical correlation (57) 
which is claimed to hold for various systems covering a wide range 
of physical properties of liquid, 

The measured bubble sizes are smaller than would be predicted 
from empirical estimating procedures in the literature but these 
were developed only for the air-water system. Yoshida and Akita 
(61) noted smaller bubble sizes in electrolytes and this could 
possibly be expected with the organic system dealt with here. The 
interfacial areas predicted from Calderbank's equation for interfacial 
area in agitated vessels, are smaller than the experimental values 
obtained in this work, though for the lower flow rate of oxygen (i.e. 
85 ccs/sec), the difference is small. Hence one must be careful in 


using the available correlations for systems other than the air- 
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The best-fit values of physical mass transfer coefficient Cee 
obtained from the model predictions and the experimental data for the 
acetaldehyde concentration profile in the sparged reactor, seem to be 
in reasonably good agreement with those predicted from Calderband's 
equation (5.14) for bubble swarms of average bubble diameter less 
than 0.25 cm and from Hughmark's equation (5.9) for bubble swarms. 

For the ideal operating conditions, i.e. well-mixed liquid and 
pure gas, used in this work, a simple model suffices and can be used 
to predict the performance of the laboratory column with a good degree 
of accuracy. However, for an industrial column air would probably be 
used and a significant depletion of oxygen would occur. In such a 
case, the spatial variation of bubbles over the reactor length and 
the decrease in partial pressure of the reacting gas component must be 
considered. Equation (4.18) then provides the gas-phase model for a 


semi-flow reactor. 


CH : 
Kaa, RT aot) "TY) de —aOl oe a oped > CRB) 


and can be solved for Yo(2), mole-fraction of component C in the gas 
phase, as a function of reactor height z, and then average values of 
ans vs and Mie must be used in the liquid-phase model of interest to 


solve for c(t) or b(t), as mentioned in chapter IV. 
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LIST OF SYMBOLS 


Ay» ans a = empirical constants in eqn. (5.7) 


a = ratio of surface area to volume of a single bubble, te 
a, = average value of a, 

a = interfacial area per unit volume of Neco bie ol ee 

a' = interfacial area per unit volume of liquid-gas mixture, L 


A = integration constant in eqn. (3.18) 


as by = empirical constants used in equation (5.9) 


b = concentration of liquid-phase reactant, mL ° 


Do = Dulk-iiquid value of b:,; mL © 


b = volume % of acetaldehyde in liquid 


mg = volume % of acetaldehyde in the liquid feed. 


: ; ‘ ake ao 
€ =.concentration of reacting gas component in the liquid, mL 


Co = gm-moles of 0. per cc of solution, mL ° 
2 
: : : ; pontine a 
Cc. = Liler Gencenuration of. C an. thes laqguid,, mL 
C' = equilibrium value of C, mi ° 


F -3 
ae = interfacial value of C, mL 


-3 
bulk-liquid value of C, mL 


a = 
O 
os = volume % of 0, in the gas-mixture. 
Zz 


qd. = bubble diameter, L 

d = bubble diameter at the orifice, L. 

qd. = bubble diameter of i number of bubbles in a sample, L 
i 


bubble diameter read by microscope (Uncorrected), L 


ay 
iH) 


diameter of the orifice, L 


Q. 
il 


d = volume surface mean diameter of bubbles, L 
VS 
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tower (column) diameter, L 
diffusivity. ee 


diffusivity of the reacting gas component into the Ligue, Lt 


2 etl. 


diffusivity of the liquid=phase reactant, LT 


activation energy 

fugacity of pure gas at initial conditions taken as 
1 atmosphere. 

fugacity of hypothetical liquid at 1 atmosphere. 
a variable defined by equation (4.15) 


acceleration due fo gravity, ijt 


volumetric flow of gas, Te 


molar gas flow rate at inlet, nie 


molar gas flow rate at outlet, ar 
Henry's constant 

a parameter used in the integration of the rate expression, 
overall gas-phase mass transfer. 

first or pseudo first order reaction rate constant, a 
reaction rate constant, units depend on the order of 
reaction. 

frequency factor 
i 


liquid-side mass transfer coefficient, ie 
gas-side mass transfer coefficient, moles/area-time-pressure 
reactor length, L 


length of a bubble, L 


latent heat of solvent, cal/gm 


i. 
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= latent heat of solute, cal/gm 
molecular weight, m 
= molecular weight of solvent, m 


instantaneous rate of chemical absorption, Fy ies le 


average rate of chemical absorption, aie We 

= total number of bubbles per unit volume of liquid 
order of reaction 

= number of bubbles of diameter dis in a sample 
partial pressure, force/area 

total pressure, force/area 

= (P-V,)5 force/area 
= critical pressure, force/area 

= agitator shaft power during dispersion. 
stoichiometric coefficient 
volumetric liquid flow rate, ee 
kinetic rate of reaction 

Universal gas-law constant 
Dackwert's model parameter, T 
= ae Qonst all k) = a 

time, T 
= total time elapsed from the moment the surface element 
considered has been brought to the surface, T 

= peaction time, T 


=nansrusiom times T 


O 
temperature, K 
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TY = reduced temperature 

T. = critical temperature 

UL, = characteristic of the liquid velocity, tre 

Ugg superficial gas velocity, fy te 

Uo. = Superficial liquid velocity, ee 

U, = rise velocity of a bubble, ur + 

U, = slip velocity, LT” 

OT AG = linear liquid or gas velocity, bas 

i = terminal velocity, ur + 

Ne = volume of agitated mass, ie 

VE = molar volume of the solute, he 

wey = vapor pressure of AcH at reactor temperature, force/area 
Vv, = vapor pressure of AcOH at reactor temperature, force/area 
Vi, = volume of liquid in the reactor, if 

ihre = volume of a bubble at the orifice, a 

ue = volume of a single bubble in the reactor, i 

W = a constant defined by eqn. (4.2) 


w = width of a bubble, L 

X = mole fraction in the liquid 

x = average inside diameter of the column, calculated 
photographically 

x' = solubility of gas in liquid, g-moles of gas/g-mole of liquid 

Y. = mole fraction of component C in the inlet gas 


Z% = distance along the axis of the reactor, L. 
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Greek Symbols: 


loa 


the aerated liquid height 


the sclear liquid height 


a variable used in equation (3.20) 
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pore diameter, L 


65 & 8, = solubility parameters 


€ = fractional gas hold-up without liquid flow 
e' = fractional gas hold-up with liquid flow 
a= (Ges) )/UC8 Cu) 
° O 

U =viseesity , mutt 

aoe Se 
Re Vuseosity of continuous phase, ML “T 
ce kinematic viscosity of liquid, Barts 
Py = specific gravity of liquid and gas, respectively 

> 
ee specific gravity of dispersed phase 


o = surface tension of liquid, force/length 


¢, = association parameter 


Dimensionless Numbers: 


ies 


=, 
i 


Froude number 


i 


Nae = Peclet number = dU, /D 


Nee = Reynold's number = DUp/u 

Nog = Schmidt number = u/pD 
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Nop = Sherwood number = Ky, d/D 
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APPENDIX A 


DERIVATION OF THE CHEMICAL REACTION RATE EXPRESSION 


Assuming that the reaction in a stirred-tank reactor can be run 
under the kinetic regime, it is desired to derive the rate expression 
for a batch gas and liquid reaction: 

CH ,CHO+0,, > CH COOH 
The pressure in the reactor is maintained constant and the 


reaction is assumed to be a first order in oxygen and also a first 


order in acetaldehyde concentration. 


Hence, dD oH ee 
ae | etaoR AcH*Po, 


It is assumed that the reactor contains only oxygen, acetaldehyde and 


acetic acid. 
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where X. = mole fraction of i in the liquid 
i 


Va = vapor pressure of AcH 
V, = vapor pressure of AcOH 
Dy y= volume % of AcH and AcOH in the liquid mixture of 
e 
AcH and AcOH. 
i 1G@0-b. 4) k.05 

D poy s+ 78) D aon! 78) ( AcH 
Lise SB aan uk ME np 
AcH 44 


(specific gravity of AcH = 0.78 and specific gravity of Aco = 1.05) 


It is assumed that there is no change in volume when acetaldehyde 


and acetic acid are mixed. 
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Substituting for Po in the rate equation, one has 
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AcH _ -2 
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= 40°C, V, = 26 psia, V, = 0.657 psia 


At-a constant temperature, say T 
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Integrating the above equation, one has, 
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In general at any temperature, 
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Determination of Solubility of Oxygen im Acetic Acid :— 


Solubility of oxygen in acetic acid is not available in the 
literature and the following equation proposed by Shair and Prausnitz (73) 
was used to estimate it, using the solubility of nitrogen in acetic acid 
(to get the solubility parameter for acetic acid) at 25°C and 


aLmesphere pressure as 0.109 ce Ni per cde liguid (75): 
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tc 2 2 
ae 7 ae exp { RT 
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where £5 = fugacity of hypothetical liquid at 1 atm. 
- = fugacity of pure gas at initial conditions taken as 1 atmosphere 


P = critical pressure 


53 65> solubility parameters 


e, = association parameter 
le aes molar volume of solute 
x, = solubility of gas in liquid. 


Acetic Acid 


CH. COOH = 05 gms/ecs. M =.60<05 
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Nitrogen 
cen, 
Solubility = 0.109 -——, |. 
ec liguid 


- Solubility of nitrogen in mole fraction, 


0.109x60.05 - 
x's SS ee on = 2.79210 
PIES PUT NOR OP Rca 
Ree A7 208.15 =1599 5 ca l/anenele 
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where 65 is the solubility parameter for acetic acid. 


For Oxygen at Doae 
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ie) 
For 0, at 40°C , RT =y1.987x313 = 622 
Lie ire 
we 5 orang gles Aa O2's f, = 5.8x49,7 = 288 atm. 
iL 2 
= 635, 74) 
= 288 eee 
Xe exp { Shak 


en 
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Henry's Law: P = xO 
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g mole AcOH 


g mole 0, 


g mole AcOH 


At 34.7 psia 
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Solubility of oxygen at 34.7 psia and 4o°c 
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Diffusivity of Oxygen in Acetic Acid:- 


The following equation, proposed by Sitaraman, Ibrahim and Kuloor 


(74) appears to be the best for estimating the diffusivity: 
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where D = diffusion coefficient, on sec 


HU, = viscosity of solvent, cp 

= molecular volume of solute, ccs/g-mole 

Mo = molecular weight of solvent 

L, = latent heat of vaporisation of solvent at normel boiling 


point, cal/gm. 
Leg = latent heat of solute, cal/gm. 


Here,T = 273+40°C = 313°K 


u = 0.98 cp (41) 


Le = 96.75 cals/gm (41) 
Vu = molecular volume of 0, = 2 em? /gm 
L = 50.8 cals/gm 
SO . 
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Dy ~ AcOH eM ectop al cm / 


Determination of the Reaction Regime 


T = 40°C, P = 34.7 psia, Catalyst concentration = 140 ppm 


K. = 3.46x10.° 1/psia.min 
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Initial concentration of acetaldehyde = 10 vol % = 0.00177 g mole/ 


eq solutions 


H = 1.423x210° psi/gm-mole 0,/cc Solution. 
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*, Initial rate (pseudo first order in 0,) 
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(1) For Ke = 0.0065 “em/sec, TD cs (0.0065) 2 = 0.554 sec. 
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(2)) For og = 0.0433 om/sec, th = To70n33) Qo “sec =O 70125 sec 
4 th << tp 


Hence, the operation is in slow reaction regime. 
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Rk PROGRAM =124 
VATN 
MAIN 
RUNGE KUTTA METHOD 


DIMENSION DES (20)s VAR (11093) 
COMMON FY{(10) s¥(10) s¥Y¥(1092000) sXX(2000)sD(1054) 5 
L1DD(5) sxXeNeIeco 


NOMENCLATURE REQUIRED TO USE THE DECK 
mee STEED SIZE 
N===THE NUMBER OF EQUATIONS TO BE SOLVED 
Mee==THE NUMBER OF STEPS TO BE TAKEN 
Mie==THE FREQUENCY OF PRINTOUT 
COm==THE VALUE OF MASS TRANSFER COEFFICIENT 
*INTERFACIAL AREA PER UNIT VOLUME OF LIQUID 
DESIRED MI=1lsPRINTS OUT VALUES AT EACH STEP, 
MI=29AT EVERY SECOND STEPsETCe 
DES(J)="=VECTOR SET ASIDE TO STORE TITLE 
OF THE PROBLEM 
VAR(J)—==MATRIX SET ASIDE TO STORE THE 
NAMES OF VARIABLES 
DATA CARDS= 
FIRST CARD=PUNCH THE TITLE OF THE PROBLEM 
SECOND CARD=PUNCH THE VALUES OF HsNoMoeMI 
(FORMAT F10059315) 
THIRD CARD===THE VALUE OF CO(FORMAT F105} 
FOURTH /CARD==<=THE NAME OF THE INDEPENDENT 
VARIABLE AND ITS INITIAL VALUE(FORMAT 3A4s 
Fl0e5) 
FIFTH CARD===THE NAME OF THE DEPENDENT VAR= 
IABLE IN THE FIRST EQUATION IN SUBROUTINE 
FUNC AND ITS INITIAL VALUE(FORMAT 3A49F1005) 
6 
6 


® 
LAST CARD<-=THE NAME OF THE DEPENDENT VAR= 
IABLE IN THE LAST EQUATION IN SUBROUTINE 
FUNC AND ITS INITIAL VALUE(FORMAT3A49F1005) 


READ (S591) (DES(U)s wv = 1920) 
FORMAT (20 A4) 
FORMAT (Fl0e5s9 315 
READ (S92) He Ns M 
READ(5s3) CO 
FORMAT( F10e5) 

NK. = N+ 1 

READ (594) (VAR(1l9K)» K = 193) X 
FORMAT (3A4»s F1005) 
FORMAT(TO's X=" oFl2e9) 

DO 5 J = 2» NK 
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R=K» PROGRAM =125 


MAIN ees (CONT'D) 
wu ae gi ae? 7 
> READ (594) (VAR (UsK)» K = Lb3)» MA WOT 

XX(1) = xX ; 

DON Z)!) me 2 

-DD(2) = 065 

DD({3) = 05 

DD{4) = le 

DDO(5) =1le 


COMMENCEMENT OF THE R=kK PROCEDURE 


DO 6 I=slaeM 
DO 7 L=loeN 
Peres YX GLa) 
DO 8 Jales4 
CALL FUNC 
DO 9 LelaeN 
9 DilLel) =H*FY(L) 
DO1O L=loN 
LO YCLIRVYY(LsI)+Dilesu)* DD{U+l]) 
8 X= XX(T)+H*DD (U+1) 
23 KeI+1 
XX(K) =X 
DO6 L=loen 
* YY¥{bLe9K) =Y¥¢too!) 
DO 6 J=ls4 
6 Y¥(LoK)=Y¥(LoK)+D( Loud) /(60*DD(u)) 
DO 16 K=le™ 
T6 YY(NsK)RVYY (Ne K)/06000177 
24 WRITE(6920) Xx(1) 
KOUNT=MI 
WRITE(6011) (DES(U)e Yvele20) 
11 FORMAT (LH19//916X92044/) 
WRITE(6026} CO 
26 FORMAT(!'O's15Xes'MASS TRe COEFF*#VOL INTERFACIAL AREAs! 
*9F1005) 
WRITE(60912) ((VAR (Lod) sJ=193) sLelsNk) 
12 FORMAT('O'»15X910(3A494X)/) 
DO13 T=loM 
TFLKOUNTsNEsMI) GO TO 13 
15 KOUNT=0 
WRITE (6914) XX( 1) 9 (YY(LoT) sLeloenN) 
14 FORMAT(15X98(F11e6894Xx)) 
13 KOUNT= KOUNT+1 
STOP 
END 
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_ R=K PROGRAM =126 
SUBROUTINE FUNC 


SUBROUTINE FUNC 

COMMON FY(10) 9¥(10) 9¥¥(10 92000) 9XxX( 2000) sD(1094) 9 
1DD(5) sX9NeTsCo 

FY(1)=CO#(14e6E= “4m 600e%#V(1) )=(2046F+3)%Y (1) "Y(2) 
FY(2)s=( 46¢92E+3)#¥(1)"Y(2) 

RETURN 

END 
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APPENDIX B 


TABLE-Be2 


THE EFFECT OF AGITATION ON REACTION RATE 


gg] 
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Pressure = 24.7 psia 


Temperature = 50°C 


ry 
i 


AcH?= 50 ecs, ACOH = 450 ccs, 


Catalyst = 17 mgs = 140 ppm (Mn** in AcH) 


Run Number Stirrer RPM Time Volume % AcOH 
(Mins ) 
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TABLE B.1 CONTINUED 


Run Number Stirrer RPM Time Volume % AcOH 
(Mins ) 
3 2910 0 JOR 
5 93.9 
LO OSip 
LS head 
Zo 97.9 
4O 98.5 
60 LOO. O 
4 600 0 90.5 
5 9.3 
LZ Je © 
ley 92.2 
25 92.4 
40 33,6 
60 ii AL 
90 WOKE O 
5 1600 0 90.5 
5 S72 
cae 98.8 
18 99.4 
25 100.0 
4.0 LOORO 


TABLE B.1 CONTINUED 


Run Number Stirrer RPM Time Volume % AcOH 

(Mins ) 

6 UTES @) SO 

(Reproducibility 

Run ) 5 94,2 

ING 96.6 

IL 98.0 

OS IOO} LG, 
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TABLE B.2 


THE EFFECT OF AGITATION ON REACTION RATE 


ae) 
i) 
WwW 
4S 
— 
4G) 
72) 
P- 
re¥) 


R= 30.-C 
AcH = 50 ces, AcOH = 450 ccs 
Catalyst = 17 mgms = 140 ppm 
Run Number Stirrer RPM Time Volume % AcOh 
(Mins ) 
7 1600 0 S055 
5 Cheri) 
10 ene) il 
aes 99.4 
20 100.0 
8 2S 5 0 90.7 
5 98.2 
10 98.8 
Ls 99.4 
20 So e7 
25 100.0 
9 2190 @) Ono 
S) oy ell 
10 98.8 
WS heal 
20 99.4 
29 99.7 
30 100.0 
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TABLE B.3 


THE EFFECT OF AGITATION ON REACTION RATE 
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Ta a501C 
AcH = 50 ccs, AcOH = 450 ccs. 
Catalyst = 17 mgms = 140 ppm 
Run Number Stirrer RPM Time Volume % AcOH 
(Mins ) 
10 21.35 0 90.5 
3 98.0 
5 99:64 
10 99.7 
iS 100.0 
Li 2910 0 90.7 
3 O06 
= 99.4 
ao) 9987 
ahs) 100.0 
12 1600 0 90-5 
5 98.2 
LO 93.8 
15 Seegc 
20 gS3.5 
Ze 100.0 
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TABLE B.4 


THE EFFECT OF SAMPLE SIZE AND AGITATION ON REACTION RATE 


‘Ud 
WN 


24.7 psia 


50°C 


rH 
tt 


Catalyst = 140 ppm 
ACOH = 900 ccs, AcH = 100 ccs 


Total Reaction Mixture = 1000 ccs 


Run Number Stirrer RPM Time Volume % AcOH 
(Mins ) 

13 1600 0 SOs 
5 heirs) 

LO Neko AL 

ILS) 98.8 

25 100.0 

40 100.0 

14 2135 0 90.5 
5 96.0 

10 96.5 

aS Chane i 

30 97.5 

45 98.6 

90 LOO .0 
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TABLE B.5 


EFFECT OF TEMPERATURE ON REACTION RATE 
P = 34.7 psia 
RPM = 1600 


Catalyst Concen. = 140 ppm. 


Run Number Temperature Time Volume % AcOH 

(Mins ) 

15 40°C 0 90.5 

5 95.0 

10 lio) 

ahs eps 

20 Siebel 

ee 99.4 

30 99.7 

SS 100.0 

16 20°C 0 90.5 

5 O10 

10 94.9 

AES 96.8 

20 iarodl 

30 Ske As) 

40 99.0 

50 99...7 

60 100.0 


HOoA # omuLoV 


TABLE B.6 


THE EFFECT OF CATALYST CONCENTRATION ON REACTION RATE 


P = 34.7 psia 


iia 40- 
RPM = 1600 
Run Number Catalyst ppm Time Volume % 
(Mins ) 
aly 83 0 90.5 
Ss 94.0 
10 96.0 
tlk) Sell 
20 Oe! 
30 98.2 
35 98.8 
40 99.4 
45 99.7 
50 ANGORA) 
Le 249 0 90.5 
5 Wiehe) 
10 99.4 
iS oe ey) 
20 100.0 
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TABLE B. 7 


CHECK OF THE REPRODUCIBILITY 


Hehe) 


REM = seo Oo 

Run Number Time Volume % AcOH Volume % AcH Comments 
(Mins ) 

13 ¢) 90.5 9.5 Check the 
Reproducibility 
5 94.8 Sy of Run Number 
WS 

10 Qed Zand 
eS 98.6 dee 
20 99.1 0.9 
25 99.4 On 6 
G5) ONO, O QO 


Ae 2a 


etnsanto? HoA # ompLoV HOOA & omuLloV 


ait ‘oso 2,8 2,08 
yiilidioubouges 
ssdmut oto Se 8 He 2 ~ 


ef 


oO 


i HOC 


P 


34.7 psia 
Catalyst = 140 ppm 
b = 20% by volume 


RPM = 1600 


TABLE. Bs 73 


THE CATALYST ACTIVITY 


SS SS SSS SS 


Run Number Time 


20 0) 


21 0) 


Volume % 


Comments 


AcOH formed above in 
Run No. 20 was used, 
Fresh AcH was used 
and make-up catalyst 
was added to bring 
the jcatalyst comen. 


to 140 ppm. 
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TABLE B.8 


THE EFFECT OF AcH AND 0, CONCENTRATIONS ON REACTION RATE 


P= 2h.) petagele=n50.C 00 See Aonsepes ae 


, 1 = 34 psia,P! = P-V, = (24.7-1.045) 
=23.655 psia, Veal: = 32.955 psia 
: =2 
. © T.O15x10° “(32.955) = 0.33% pata; b , . = 9.5% 
AcH 
P'-0,334 5, b°AcH 
Rate Equation: I1n( 5 re ee 
H _ . te 
Ac P'-0.334 DCH 
fo) 
ve _ o05 eur 
SWieeine 23.655 —- 0.334(9.5) 
AcH 
(P'-0.334 b ) 
or, in 4 5 _ CP SS) 8 mm a 
AcH 


Catalyst = 140 ppm 


22.0050 
22.987 
23%:L00 lO 2 
26.300 11,936 
23.421 15.390 
23a Didrds 27.047 


~08351 23.571 43.368 


0334] 23.621 | 108.650 


ATAA MOLTOAAR HO, GMOTRAR EMS? Mth Heh 9 
i | A ae 
(auy.t-0.88) = .V-d = 4,5ieg = ,V .pteg Gul = ot oh 
plaq ee = pe 


o2.0 + . °d gekag BOE.0 = (088.5 


HOA +) 
Hoa d d Hee. On'a ¥ 
a ee ee. = — . : 4 s 
mA £. * -' HoA 
Hoa #EE.O re 


- ft ree 

oa 2 mes _ e os : ‘ : - > r. 

aaa CES - med 6 
ie as 


{ d #€&.0- 'g) a 7 
Gogh = (3R.0) » AR nemen fil gt 
Hos 7 iP a 

mqq OL = taylesso 


£8.94 
evs iL 
S$ PBL 
2236S 
are. V0E 
$8. 426 
ey.ese | 
OL .eve | 


¥9,s62 | 
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TABLE. B.9 


TRAE EREECT OF AcH .and 0, CONCENTRATIONS ON REACTION RATE 


P = 34.7 psia, T= 50°C =P! = 34.7-1.045 = 33.655 psia 
e) 
e AcH = OO 
B AcH 3 er = 0.312 
te ence oe ete 33.655-3,17 
AcH 
P'-0,334 b 
ln { = AcH (9.03) } = KtP! 


Catalyst = 140 ppm 


ie Aine) Bacon) Aci Be il etLee t.Pt 
AcH P=, 33 

Vi 2.895 67 2o 0 
33.154 6.895 Ah, Byes 

OG DIS 8.646 Deo woeo 

88) Sos ALL, ES} SS86n100 

33,455 1 5 BSNS ISNA) os '5) AL 
Gn See! 26.145 BOW, 825 
Beha, Se 41.896 Sy eg alee hs: 

33.622] 104.490 639.44 


| a 
Pe Whe Te White 


PAR WOLTIAda WO ¢ 


—s =f : ‘ ot, : 4 
eing 22a.8e = eH0.f-T.He = Sh OO 


$ 
a 


oF 


(Sib. 


Hos 


eeB.< 


208.8 


ore: ¢2 
ane.6ar 
age.ees ava.8 | 
26.a8€ p8e.fL | 
£2 ,.eil ae WL 


258 ..uoe) aer.as 
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TABLE B. LO 


THE EFFECT OF AcH and oO, CONCENTRATIONS ON REACTION RATE 


P = 44,7 psia, T = 50°C, 140 ppm, 2135 RPM 


= 5% 
Pig € 22 

12) 

: AcH . 5 » 5 - (119 
BI. 394 H° 43.655-.334(5) 41.985 ~ 


Fy Newe 


so2e 461.5 


O.loe Wee c521 306.0 
0.0668 | 43.588 424.0 


q.t 


i(err.) 0% 
L 0 


‘BA 
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P: : as a a ey ” , 
ATAA MOITOAZA MO ako TTARTHSOHO 
Nes P é A 


eft. 
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TABLE Beil 


THE EFFECT OF TEMPERATURE ON REACTION RATE 


P = 34.7 psia, Catalyst = 140 ppm, T = 40 C, Vy = 26 psia, Vy = 0.657 *psia, 
V,-V, = 25.347 psia, P' = 34.7-0.657 = 34.043 psia 

5 =2 -2 

emo ete “Ve Vi) = Ore x 0 AC2o. cur). = ON 258 ° 

ae b 
‘ AcH = ee 
ae gers nee 31.493 
AcH 
Pr=0.258 b é 
in { = pe (d, 302) = Kee ee 
AcH 


1S S23 700 Ue lee iks: 
33.140 340.43 
Shei Morel 510.645 
33.811 680. 86 
33.889 851.0795 


33.966 


121.290 


ae eae 
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dates i (an 
io she any - ‘ ‘oly uA pare - 


»steq Yed.0 = ¥ eeteq as & Fa oe ou eT ema | 
steq s#0.#6 = yaa.0-¥ He = snai 
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eerste = Ci aise <ts 
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ais a0 'g.2% = 1 (S08 +0) ® . 
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avo, £28 ba;S e868 66: | 
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141 


TABLE. B..12 
THE EFFECT OF TEMPERATURE ON REACTION RATE 
Po= ops /,.psia, Catalyst = 140 ppm, T= a0 c. V5 = 20 psia, V, = 0.387 psia, 
VV So 6 lowpsias Fig = P-V,, = 34,7 =~ 0,987.= 34.313 psia 


ee 1.015x10°“(V,-V,) = 1,015x10 “(19.613) = 0.199 psia 


fe) b- 
b = 9.5; AcH : g.5 - 
noe P.'-0, 198 b 2S. 
AcH 
P1=0,199 b 
in ne (0.294)}= K.t.P! 


»5229 Tee .0 = 


a 
Bieq PPL.0 = (€L8s0L)> OLXELOWL = (g¥- pW) omnes. 6. 7 


o 
2,2 A ina © 2.0") de” ; 
weS-0" CE - elewe ‘ 
“em - cere ee 
| a 
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TABLE B13 
THE EFFECT OF CATALYST CONCENTRATION ON REACTION RATE 


P = 34.7 psia, T = 40° o, Catalyst = 83 ppm, V. = 26 psia, V, = 0.657 psia, 


a 2 


V,-V5 = 29.347 psia, P! =°9394.7-0.657 = 34.043 psia 


=2 x 
= « Le015xl0 Veale, = 1.015x10 2(25.347) = On208 


Bp 
AcH a.5 
O . . = === = 0.302 
bach = 9:52 *"* Bt_9. 958 4° 31.493 
AcH 
'- SE 
in—4 pa (0.302 )}= K t P! 
AcH 
Run Time | 5 b -258 |Py = Po (.302) Py (.302) 
No. | (Mins) os AS os ¢ + | int | ste?! 
AcH | -P'=.258 D b 
5 AcH AcH 
AcH 
17 94.0 6.0 1.548 | 32.495 1,635 0.492 170.9 

96.3 ou 0.954 | 33.089 2.7100 0992 siame 
elisa Pe 0.748 4 32.295 3.467 1.24 hee 
ly 2.3 0.593 | 33.450 4,391 1.48 ae 
982 1.8 | 0.464 | 33.579 5.633 tere nee 
98.6 ie | O.361 | 033.682 7.265 1.98 1699 
98.9 elev 0s 2820) ecceuie 9.268 2.22 1160 
90.654) 0.05 |) (Oiled 4) “one 75 15.738 Vile 


} 1365 


,bteq $23.0 


(SOE). ot | (508.9 ot] = “a | 
bi {ie er ae sas5- | 
HoA rat | 
2. ORE S2u,.0 260.1 au SE } 
a ee | 
aie Wolk 9 |) Naw.e aes, bef 
ea aii res.e | oew.ee | 
Pore ee a eva.ce | 
Pe ee 2as.y | SBase€ | 4 
ret BES | PS pase =| av ee | 
sank a jars } sere | ateuet. 


an ty sieq aS > = wv - neg 


eteq eno.ve = heey aul 
Bes.0°= (THE, pia 


“a y 
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TABLE B.14 
THE EFFECT OF CATALYST CONCENTRATION ON REACTION RATE 


P = 34.7 psia, T = 40°C, Catalyst = 249 ppm, P' = 34.043 psia 


5° 
fe) 7 AcH O45 
b = 25 Gian = —— = 0.302 
AcH fe) 315493 : 
P'-0.258 r 
: AcH 
R : ; p ss p p 
un Time b OOH DoH 0.258 0, 0,,(.302) 0,,(- 302) 
No. (Mins ) b ee in—————— t .P! 
AcH Pbi—.258 b b 
b AcH AcH 
AcH 
18 @) 90.5 9.5 25 Shik, Sy ehs} ab ¢) 0) 
Did 915.5 4,5 dhe Aus) B25 Site! 22D 0.79 SHS) a Z 
5 9725 DES Oey sis 5 Seis 4.07 ee SIS ALGO, 
Hos 98.5 aes} OG sks Sha (5 NE 5s nba sche PRY pO) 
10 99.15 Onc Og 22 chele V2.6! eS 2S 340.5 


ats) 99255 0.35 @..090, 733.953 29.4 3.38 307,10 


c OE 


0.0L? 


pee). Le 
ese.sé ai.f  &.# eee 


, | | 
Bee.ce e400 gat ee * | 
ged .GE V8E.6 if 2.88 


~ 
ww 
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4) 
Se) 
i) 
3 
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ae! 
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TABLE B.15 


COMPARISON OF THE PREDICTED AND EXPERIMENTAL RESULTS 


P' = 34.043 psia 
Kp = Seno Veo Ge | ah pete. uu eeen 
P'.K. = 0.118 


R 


Calculated west 


Po (302) 
Experi-= Z 
Source ln ae SE ES 
mental DoH 
ryt — (Mins ) 


(Mins ) 


iota, ING. Is 


and 


table B.11 
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Obser- 0 
vation 
No. 


TABLE B.16 


OXYGEN FLOW RATE VS FRACTIONAL GAS HOLD-UP 


20 psig, T = 40°C 


5 Flow 


Rate at 

40°C and 

34.7 psia, 
ecs/sec. 


207.0 


CM 


= f  G 
cm cm 
LOZ EO. yor 
103.40 1. 800 
1062/81) Sa Si 
110.39 8.800 
tS. 60: 12,000 


0.0174 


0.0486 


0.0796 


0.1058 


0.247 


0. 343 


0. 806 
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TABLE. e/ 


COMPARISON OF THE EXPERIMENTAL e¢ WITH THE VALUE GIVEN BY HUGHMARK'S EQUATION 


e 9) 
Observation U ae d2yi/3 
SG Pr, oO 


No. 


) e(Expt) 


(from graph) 


e' (Hughmark ) 
(from Hughmark's 
graph) 


e (Hughmark ) 
(back- 
calculated) 


PEO .0 


a620.0 


( 0L.0 
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TABLE Be 


INDIVIDUAL BUBBLE DIAMETERS 


0, flow rate = 334 ccs/sec, magnification = 10, column diameter in DNOtG = 
[Seg eee 
os COMP. ac70 
1 Z 6) 
a a a 

18 iS by BROS) 1.94 15.04 PAS! AG 
£2 10 ae AS LK) 45.63 ic iyce 
12 se uy 22 dace 7.60 8. 74 
0 10 38 10 ile iS) DOs. L6 D7 268 
14 ad 1 ez ais 1.49 PIE 3.30 
14 10 4 LZ eS 1.44 2007 Paposhs| 
7 7 a 7 0.805 0.64 Oey 
ales) 13 2 14.3 1.64 5. 36 8.79 
13 13 ak is 1.49 eiZece 3.30 
20 20 4 20 eae) 21.16 48.66 
10 8 7 Oven? 7.98 8.53 
20 8 1 14.7 1269 acs 4. 81 
5 10 9g LEC ial Liou 20S2 30.98 
45 als 24 Ls Lai2s 70. 80 ibe Wa) 
10 9 a Gy GHe Seadit 1,23 236 
20 15 7 18.2 2.09 30.52 63,78 
20 12 3 16. 86 1.94 ii, 28 feces 
LS 12 ie 13.92 1.60 2456 4. OD 
20 18 2 “19.3 2.22 9. 84 21.84 
18 8 a Load OS SHTAS. 2.46 3. 86 


G Cont dy 
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TABLE B.18 CONTINUED 


Z 3! 
1 W Number of 4 aides Ti. 
Bubbles, n, “b “b i% ip 
20 7 tL NSS I ats Zeon 4.89 
20 10 2 15.9 “les (ks 6.68 Re ee 
30 igs it DXEVA I} 2/4 Teta) PANG SS) 
30 20 alg 262 8, (02 Orlez 27,54 
doh Ole COLL y=570). 46 
29 7/0..46 Senate 


vs 348.31 
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TABLE: Be 19 
INDIVIDUAL BUBBLE DIAMETERS 


0, flow rate = 215 ccs/sec, magnification = 10, Column diameter in Puota = 2, 


Bio ors 

a Lo one 

1 a Number of qd d na n.d* 

Bubbles, n, 

‘a 5 34 5 784 20.40 Me am 

6 6 ay 6 - 94 0.88 Oe 

8 de: 8 6.3 heh) 7.68 7.60 

8 5 a 628 APP Olezs) eee st 

8 8 8 8 aR ks) 12.48 15.00 
10 8 ak 9.3 1.46 Zeid 3.10 
ANS) a 40 8.6 L530 90.0 135.00 
LO 10 2 10 ows 4.92 TZ 
Le 8 ail Oye 1.60 Poo 4,09 
LO 10 alae 11.3 hel] 34.43 60.94 
20 4 10 daly? A tek 33.40 Glee 
LO 6 i 8. 44 Loz 1.74 2329 
10 5 ak URE 1.24 1,538 APR o)2) 
be dee i 12 1688 3.00 Gabe 
as) 10 an ilies spl 2.05 4. 20 “shes oll 
30 10 2 20.8 3.26 21.24 69.24 
30 15 a 23.3 3.68 13.54 49,82 


E=255.78 2=451. 66 


Ce a ee er 


451.66 
= SSE Se Hagley 
esse 


"“¢ i= oFedg ar poses nme £03 ger = ais 


[8 ,2£ 04 ,0¢ 
$2.0 68,0 
7m) 23,) 
Pare See 
ya el Biles 
of EL ,8 
OOVeEL 0,08 
$3.7 oni 
og. ea 4 
#2,08 Evie 
Si<fa@' ° OviEB 
es.s uv ir 
ea. bat 
$8\3 S8x& 
{3.8 OS 44! 
Se I 


— a 


a 7 ca one ae 


ae 


vay. c 
FE 4 ra) 
26e, rs 


CS ak = 
Slat €.e 
Oc.L is 
Te. OL 


Yat EEE 
Be.f VLE of 
SE,I -#H,8 t | 
WS.f .SB,Y £ 
68.1 Sr . ) ie : 
20.8 f.€L iy 
OS.£ 8,05 s 
‘ ‘ 
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TABLE B. 20 
INDIVIDUAL BUBBLE DIAMETERS 


0, flow rate = 177 ccs/sec, magnification = 10, Column diameter in photo = 2.50", 


' 2 3 
ee eh a agen Te ae 
8 8 2 8 is LOY 2.436 2,682 
8 6 a 7226 ee O) L002 £008 

LO 8 5 SNe) e283 8.23 Ls OS 
2 6 ie 9,5 de Salish Le LS Zee 
LO LO 8 10 agers: Los232 21020 
12 8 1 LOZ 1.407 MESSE) 2.784 
L2 LO 5 ins: UPI | Beare Roi) 18.941 
£2 LZ 3 lig O06 85226 13.622 
GS) 10 LO US ie EA SOW BAe Siew shshs! 
20 LO all IS. 86 Zin BB 4,788 10.476 
18 18 il 18 2.484 6.1/0 E55 3.26 
20 tis th 18.2 POS oelaih 6.305 LS Sail 


2=100. 886 2=173.494 


eee SSS a I OTT 


_ 173.494 


= es eee 7 9 cC 
vs 100. 886 Eee 


Leh. el 


3p ET i=3 


a88.00f=3 
j a 

al bed aa 7 iF a> 
ie 


oe 


~vei : , ' 

: ; y 
7 nee. ; 7 
oo \ O81 
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TABLE B.2l 
INDIVIDUAL BUBBLE DIAMETERS 


0, flow rate = 85 ccs/sec, magnification = 10, Column diameter in DINSTO: = a aO ae 


af oe 3.45 
lan as0 
2 3 
ee ae 
6 6 g 6 1.035 9.54 9,82 
8 8 47 8 Tyee 89.30 103.28 
9 g 8 g ies 19.20 29.76 
10 8 2 9.28 1.60 5.12 8.19 
10 10 31 10 tae 91.45 157.29 
12 10 4 11,28 1,945" | 16.20% yee 
13 10 2 PSG 2.05 8.40 1722 
12 il 2 12 Lagu 8.56 siiew gs 
15 10 ai 13.08 2.26 56. 10 126.78 
15 1D 3 ale ele Pag 17.20 41.28 
16 12 1 14.52 2.50 6.25 15.62 
ih ce “l 14.70 2.53 6.40 16.19 
18 10 ib 14, 80 2.55 6.50 16.57 
16 15 a 15.64 2270 7.29 19.68 
18 12 if 1SB10 pace i 7.34 19.89 
15 15 0 15.00 2.58 6.65 17.15 
20 10 3 15.86 O18 22.35 61.01 
16 14 1 15.3 2.64 6.96 ce a7 
15 15 il AC 2.59 6.70 17.35 
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pe ee PY 

7 Pua TaMALG ogtoave 
SAM MALE ote 
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TABLE, Baad CONTINUED 


1 : Number of 2 3 


Bubbles, n. a a rags aia 


14 10 lb EWP? 2.16 4.66 10.06 
20 lhe al 16. 86 2 ou 8.46 24.61 
20 LS ae sre US) es) Says 30.64 
25 LO ne 18.4 3.18 10.11 O2aL4 
20 ay) i Le.9 3.26 10.62 34.62 
244 aS) a: 20.4 eiice 12. 09 43.61 
30 he a 2202 3.83 14,66 okeeee re) 
40 12 dl: 2067 4.60 21.16 97.33 


2=488.11 2=1091.24 


een A AS A A AS EA A RE SAE PAGS SEES SEN IST ESS EEA SE 


S097. 25, 
Sa) STEED cul Le 


ag, #2 


aL LS 


we ee 
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TABLE Bs 22 
INDIVIDUAL BUBBLE DIAMETERS 


0, flow rate = 47 ccs/sec, magnification = 14, Column diameter in photo=1.70 " 


ep! 3.45 


of cis “Lo 


' 2 3 
as eee ee, el ee ee 
5 5 i 5 Oe t2 0.518 0.372 
6 6 i 6 O. 864 Seal’ Pwo 
8 6 uh 7.26 1.045 4,366 4.562 
8 8 3 8 1, 152 ele cjicul 4.586 
10 6 10 8.44 de 2L5. (14. 76 17.933 
10 8 8 9.3 1.339 14.336 ESR 
10 LO 9 10 1.44 13.697 26.866 
14 5 nh 10,04 1.445 2.088 Se OL7 
16 ) z 10. 84 1.56 2.433 35 195 
LS 6 ra 11.4 1.641 2.092 4417 
42 10 alt SAIS 2 LSPA: 2.647 4, 306 
ee Le a 2 A ERs! 22089 Da Loe 
EZ 8 At 10.2 1.468 Zoe 3. L63 
18 8 Alt 13.7 Lay? 3.888 7.667 
iS 13 aL 14.3 20059 4, 29 8.833 


2=84.969 Y=118. 381 


Se a ET 


118.381 
eee z 


oe “Sn. 969 Loos mM « 


Pesci pguiner aaa: 


.c t 
Lie.e SSs.¢ -#08).0 a v 
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#br er GEO HE LEE.L 6.2 a 
gealas = ta bE Ht of of e 
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TABLE B.23 


CONVERSION OF AcH to AcOH IN SPARGER REACTOR 
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20 psig, Catalyst concentration = 140 ppm, Liquid height = 40", 


it 


Diameter of column = 4", Orifice diameter = 0.2 om., Liquid recycle rate 


1.83 U.S. G.P.M., 0, flow rate = 334 ccs/sec; pa = 10% vol. at 40°C and 
34.7 psia. 
Ces b CoH D py volume % ) 
0 90.0 10.0 
6 a2 56 8.0 
10 93.3 6.7 
20 96.1 3.9 
30 98.1 19 
4.0 99,3 0.7 
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TABLE B.24 
CONVERSION OF AcH to AcOH IN SPARGER REACTOR 


0, flow rate = 85 ccs/sec, T = 4o°C, P = 20 psig, Catalyst concentration 


- ° 2 
= 140 ppm, 5b ee 10% volume. 
Observation Time b ” is ’ 
Wimber Gin) AcOH (volume %) AcH (volume %) 
1 0 90.0 10.0 
~ 10 91.6 8. 4 
: 20 93.0 7.0 
a 30 94.3 a5 
8 40 96.1 3.9 
6 50 96.8 Ome 
7 60 98.4 16 
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TABLE B.25 156 


ee ees i £0. 00000590 (4268705166 - 


37.9944000 0.00000617 

39.99391174 0.00000€45 4.14634228 
41599.342346 0200000676 _. 3288530922 
43.992935138 0.00000709 3262930012 
45299244690 0.00000745 3237870884 
47.99195862 0.C0000783 3213397026 
49299147034 0200000825 22 89553642 
51.99098206 0.00000870 2- 66390033 
53.99049377 0-00000918 Hr 48989236 
55.99000549 C.00000970 222231 71656 
S75 989ST721 0.00001026 2201524258 

9.589 20000 3757 
61.98854065 0.00001149 1.62719917 
63.98805237 0.00001217 1244828892 
65558156409 C.00001289 _1.28027058 
67.$8707581 0.00001 364 1. 12362862 
69298658752 0.00001442 0.97874540 
7 

73298561096 0.00001604 0.72532558 
75298512268 C.00001686 0.61694050 
17298463440 0200001766 0-52056825 
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TORES HF-585 


MASS TRe COEFF*VOL INTERFACIAL AREA= 0.01100 
TIME{MINS) CONCN-02 CONCN.-ACH 
1299996379 0.00000327 9.71738815 
$32999995327 0.00000337 9240493488 
ae99997a 15 0-00000347 9209396648 
LIN93892 23 0200000357 8278456020 
9.999861 72 0.00000 368 8247680569 
13.59980068 0.00000392 7 «86663 342 
15.99977016 0200000405 7256442547 
17.999282 84 0290009418 tnd 28429653 
19.$9879456 0.00000433 6 96637440 
21.59830627 0209000448 6.6707SC67T 
a9 0 8 
25999132971 0.00000482 6.08728409 
27.99684143 0.00000501 5. 79971218 
29299635315 9200000521 2931517200 
21.99586487 0.00000543 5223389149 
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TABLE B.26 id? 


CONCN-TIME PROFIL FE-OXIDN OF ACH IN BUBBIE REACTOR 


MASS TR. COEFF*VOL INTERFACIAL AREA= 0.03000 


TIME(MINS) CONCN-D2 CONCN.-ACH 
0.00000010 0.90000000 10,.00000000 


1.99998379 0.00000742 9. 56505699 
Bes 99SS2 7 0.00000781 8268914986 
enaawees 0.00000823 8.03003502 
1099989223 0.00000868 7238908291 
9.99986172 0. 00000917 62176775646 


11.99983120 0.90000970 6.16762447 


13.999 80068 0.00001027 5259034252 
15.99977016 0.00001089 5.03763008 
ATs 99928284 0.00001155 4251126957 
19.99879456 0.00001225 4.01299381 
21.9983 0627 0.00001299 3.54451 847 

23.99781799 0.00001377 eA AON ite 2 es Se ee 
25299732971 0.00001458 2.70293427 
27.99684143 0.00001541 2235217430 
29 299635315 0.00001625 14299571037 
31.99586487 0.00001710 1.69362545 
3309953 1659 0.90001793 1242545605 

ae ene OOOO T RTS. a PID PO TAR 
37.99440002 0.00001948 0.98619944 
39.99391174 0.00002018 Use LiaussZ 
41.99342346 0.00002082 0266302514 
43499293513 0.00002140 0.53862447 
45.99244690 0.00002190 0.43526036 

ee Pee i? 0, NON 2234. CO ASOT RIG 
49299147034 0.00002270 0.-28052193 
51.99098206 0.00002301 0.2.22401375 
53.99049377 0.00002326 0.17839134 
55.99000549 0.90002348 0.14173377 
§7.98951 721 0.00002365 0.11239636 

ee EG 2 Po OB B94 RG ee 
61.98854065 0.00002390 9.07038200 
63.9880523T C.00002399 0.05560732 
65.98756409 6.00002406 0.04389931 
67.98707581 0.00002412 0.03463450 
69.98658752 0.00002416 0.02731119 

ee) SEAS 6 OU 247 ee Oe ee ee 
73. 98561096 0.00002423 0.01696357 
75.98512268 0.90002425 0.01336366 
717298463440 0.00002427 0.01052557 
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